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Introduction 


The  work  for  this  research  contract  on  atmospheric-pressure  gas 
lasers  had  as  its  main  objectives  the  generation  of  short  laser  pulses 
through  mode  locking  and  cavity  dumping,  the  study  of  the  nonlinear 
amplification  processes  in  a  transversely  excited  atmospheric  pres¬ 
sure  laser,  and  development  of  computer  programs  to  check  the 
observed  amplifier  responses  against  models  of  the  relaxation  proces¬ 
ses  of  the  CC>2  molecule  in  a  laser  discharge. 

During  the  progress  of  this  work,  it  was  realized  that  a  simple 
yet  reliable  model  of  the  discharge  excitation  of  the  laser  medium  was 
needed  if  quantitative  predictions  of  laser  performance  were  to  be 
made  without  excessive  computational  effort  and  hence  a  certain  loss 
in  understanding.  With  this  in  mind  an  analytic  model  was  developed 
for  the  interaction  between  the  electrons  in  a  molecular  laser  discharge 
with  the  active  molecules. 

One  of  the  major  bottlenecks  in  the  generation  of  short  laser  pulses 
is  the  mode-locking  process.  Shorter  pulses  in  general  are  obtained 
through  the  use  of  saturable  absorber  mode  locking  as  compared  with 
forced  mode  locking.  The  saturable  absorber  mode-locking  process 
is  less  well  understood  than  the  forced  mode-locking  process.  In  an 
effort  to  gain  a  better  understanding  of  saturable  absorber  mode 
locking  a  new  theory  of  mode  locking  was  developed.  The  first  closed- 
form  solution  for  a  pulse  produced  by  saturable  absorber  mode  locking 
was  obtained. 


The  ensuing  sections  of  this  report  are  addressed  to  one  topic  at 
a  time  and  describe  briefly  the  main  thrust  of  the  investigations. 

Details  will  be  given  in  five  appendices.  Appendices  I,  II,  and  III 
contain  excerpts  from  the  Quarterly  Progress  Reports  of  the 
Research  Laboratory  of  Electronics.  Appendices  IV  and  V  are 
preprints  of  a  paper  that  has  been  accepted  for  publication  and  of  a 
paper  that  is  in  preparation. 

Mode  Locking  and  Cavity  Dumping 

In  order  to  produce  high -intensity  pulses  to  be  used  in  nonlinear 
amplifier  studies  in  a  CC>2  TEA  amplifier,  a  single  pulse  was  selected 
from  the  laser  oscillator  by  means  of  mode  locking  and  cavity  dumping. 
The  cavity  dumping  scheme  shuttled  the  mode -locked  pulse  out  of  the 
cavity  at  a  time  when  the  pulse  had  reached  maximum  intensity.  In 
this  way  the  full  energy  of  the  pulse  inside  the  cavity  was  obtained  in 
a  single  emitted  pulse. 

This  scheme  was  made  to  work  to  produce  single  pulses  4  ns  wide 
and  with  a  peak  power  of  the  order  of  a  few  kilowatts.  Unfortunately 
the  power  was  not  as  high  as  we  originally  hoped  because  the  losses 
of  the  mode -locking  crystal  and  the  electro-optic  switch  used  for  the 
cavity  dumping  were  high  enough  to  depress  the  laser  power  greatly 
below  the  value  that  would  have  been  obtained  without  the  insertion  of 
these  two  elements.  We  then  decided  to  select  one  or  more  pulses 
for  the  amplifier  study  from  the  mode -locked  pulse  train  emerging 


from  the  output  mirror.  In  this  scheme,  which  is  described  in  greater 

f 

detail  in  Appendix  I,  500  kW  pulses  with  a  pulse  duration  of  2  ns  were 
obtained. 


Measurement  of  Short-Term  Recovery  of  C02  Population  Inversion 

Sequences  of  pulses  selected  in  the  manner  described  in  the  pre¬ 
ceding  section  were  used  to  determine  the  population  inversion  (gain) 
recovery  via  V-V  relaxation  in  a  CC>2  TEA  amplifying  system.  The 
amplifier  was  filled  with  a  mixture  of  COz  N2  and  helium  at  200  Torr 
and  the  partial  pressure  of  CC>2  was  varied  in  the  experiment  from 
5  Torr  to  40  Torr.  The  individual  input  pulse  lengths  were  2  ns, 
the  spacing  between  the  pulses,  12  ns  The  first  laser  pulse  was 
made  to  enter  the  amplifier  ~30  ps  after  the  application  of  the  dis¬ 
charge  current  pulse.  The  recovery  of  the  gain  of  the  secona  pulse 
as  a  function  of  partial  pressure  of  C02  was  determined.  The  total 
pressure  in  the  amplifier  was  sufficient  for  appreciable  equilibration 
of  the  rotational  population  distribution  within  the  duration  of  one 
pulse.  Furthermore,  since  the  rotational  relaxation  is  mainly  a 
function  of  the  total  pressure,  and  not  of  the  partial  pressure  of  C02, 
variation  of  the  partial  pressure  did  not  affect  the  recovery  time  of 
the  rotational  population  distribution.  Therefore  the  increased 
recovery  as  a  function  of  increasing  pressure  is  attributable  to 
V-V  relaxation  processes.  The  level  which  has  an  appreciable  popu¬ 
lation  after  equilibration  of  the  vibrational  temperatures  and  which 
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can  feed  the  upper  laser  level  appreciably  within  12  ns  at  the  partial 
pressure  used  is  the  Oil  level.  Hence  this  experiment  ascertains  the 
rate  of  relaxation  of  the  Oil  level  and  the  degree  of  population  recovery 
by  means  of  this  relaxation  process.  Figure  1  shows  the  result  of  the 
experiment.  At  the  very  low  partial  pressures  when  the  relaxation 
mechanism  is  not  operative  within  12  ns  the  gain  of  the  second  pulse 
compared  with  that  of  the  first  one  has  been  reduced  by  20%  because 
of  the  population  depletion  of  the  first  pulse.  The  experimentally 
observed  recovery  as  a  function  of  partial  pressure  is  shown.  Two 
theoretical  curves  computed  on  the  assumption  that  the  recovery  is  due 
entirely  to  the  feeding  of  the  upper  laser  level  by  the  Oil  level  popula¬ 
tion  are  shown  as  dashed  curves.  The  two  curves  were  computed  for  two 
values  of  population  inversion  changes  consistent  with  the  gain  changes 
caused  by  the  first  pulse.  This  shows  that  the  major  portion  of  the 
recovery  is  explainable  by  this  mechanism.  The  fact  that  the  experi¬ 
mentally  observed  recovery  is  larger  than  the  one  predicted  on  the 
basis  of  the  Oil  level  relaxation  is  attributable  to  the  more  effective 
relaxation  of  the  lower  laser  level  with  increasing  partial  pressure,  an 
effect  not  taken  into  account  in  the  theory  which  assumed  that  even  at 
the  low  partial  pressures  the  lower  laser  level  relaxes  completely 
within  the  12 -ns  interval. 

Details  of  the  experiment  are  given  in  Appendix  II.  The  theoretical 
computations  are  contained  in  the  Master's  thesis  of  Y.  ManichaikulJ1^ 
The  conclusion  drawn  from  the  experiment  is  that  energy  storages 
available  in  vibrational  combination  modes  can  be  utilized  if  the  pulse 


to  be  amplified  is  lengthened.  The  present  experiment  indicates  the 
magnitude  of  the  effect  and  the  time  scales  at  which  such  utilization 
occurs. 

Development  of  Computer  Codes  for  Amplifier  Studies 

A  four-temperature  model  including  rotational  relaxation  for  the 
study  of  pulse  amplification  in  a  TEA  amplifier  has  been  developed  by 
A.  H.  M.  Ross  as  originally  proposed.  The  work  eventually  was  not 
supported  by  this  program  but  by  a  National  Science  Foundation  grant. 
Yet  Ross's  results  were  used,  for  example,  in  the  interpretation  of 
the  amplifier  experiment  discussed  above.  Ross's  program  is  out¬ 
lined  in  Appendix  III  and  is  available  upon  request. 

The  funds  freed  by  the  fact  that  NSF  support  was  available  to  Ross 
provided  the  basis  for  two  other  research  projects  not  originally  pro¬ 
posed.  These  are  described  below. 

Electron  Distribution  and  Lasing  Efficiency.  Computer  codes 
evaluating  the  electron  distribution  and  the  pumping  of  laser  levels 
have  been  developed  in  several  laboratories,  notably  at  United  Air¬ 
craft  Corporation.  Whereas  good  quantitative  results  can  be  obtained 
from  these  programs,  it  is  difficult  to  cull  from  them  physical  insight 
about  the  influence  of  various  parameters,  without  carrying  out  costly 
computations.  With  this  in  mind,  we  have  developed  a  simplified 
model  for  the  electron  distribution  and  the  pumping  of  the  electron 
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laser  levels  which  is  amenable  to  closed-form  solutions.  The  details 

r2l 

of  this  work  have  been  published.  Here  we  give  a  brief  summary  of 
the  salient  features. 

By  simplifying  the  energy  dependence  of  the  collision  cross  sec¬ 
tions  of  the  electrons,  it  was  possible  to  obtain  a  closed-form  solution 
for  the  electron  distribution,  a  function  of  the  degree  of  molecular 
excitation.  In  this  way  the  feedback  could  be  ascertained  of  the  changes 
in  the  molecular  population  inversion  upon  the  electron  distribution. 
Also  the  question  could  be  asked  to  what  extent  the  i-v  characteristic 
of  a  preionized  laser  (E-beam  laser)  depends  upon  the  lasing  action. 

It  was  found  that  the  i-v  characteristic  depends  only  weakly  upon 
lasing;  in  fact  when  the  elastic  collision  frequency  of  the  electrons  is 
assumed  independent  of  electron  energy,  the  i-v  characteristic  is 
entirely  independent  of  whether  lasing  does  or  does  not  occur.  The 
energy  extraction  by  induced  emission  is  offset  by  a  decrease  of  the 
energy  transfer  to  the  translational  modes  of  the  gas.  For  other 
details  we  refer  to  the  publication  itself. 

Mode -Locking  Theory 

A  major  bottleneck  in  the  production  of  short  laser  pulses  is  the 
mode-locked  oscillator.  Extensive  work  has  been  done  on  the  theory 
of  mode  locking,  but  no  simple  analysis  has  been  published  on  the 
theory  of  mode  locking  by  a  saturable  absorber.  In  an  effort  to  develop 
a  closed-form  theory  for  saturable  absorber  mode  locking,  we  arrived 


at  a  greatly  simplified  theory  of  forced  mode  locking  which  enabled  us 
to  solve  hitherto  unsolved  problems.  A  paper  has  been  submitted  to 
the  Journal  of  Quantum  Electronics  and  revised  upon  request  of  the 
reviewers.  The  revised  version  is  attached  as  Appendix  IV.  The 
r.ew  problems  that  have  been  solved  are  (a)  forced  mode  locking  of  an 
inhomogeneously  broadened  laser  medium,  (b)  forced  mode  locking  by 
square -wave  amplitude  modulation,  and  (c)  stability  analysis  of  the 
forced  mode -locked  pulse  train. 

A  closed -form  solution  has  been  obtained  for  the  mode -locked 
pulse  of  a  homogeneously  broadened  laser  mode  locked  by  a  saturable 
absorber  of  relaxation  time  that  is  short  compared  with  the  pulse 
length.  Details  are  given  in  a  preprint  which  forms  Appendix  V. 
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1.  GENERATION  AND  AMPLIFICATION  OF  HIGH-INTENSITY, 

NANOSECOND  TEA  CO.,  LASER  PULSES 
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U.  S.  Army  —  Research  Office  —  Durham  (Contract  DAHC04-72-C-0044) 

Y.  Manichaikul 

We  have  reported  previously  on  generation  of  short  nanosecond  ,>uis  •>  from  a  pin 
resistor  TEA  CO.,  laser  by  way  of  mode  locking  and  cavity  dumping.1  Single  pulses 
4  ns  wide  (full  width  at  half  maximum,  FWHM)  with  peak  power  ~1  k\\  were  produced. 
Problems  were  encountered  with  the  system  and  a  new  one  was  built.  The  problems, 
and  the  changes  that  have  been  made  in  the  new  system,  are  as  follows. 

(i)  The  peak  power  obtained  from  the  previous  system  was  too  low  for  some  exper¬ 
iments.  For  example,  to  saturate  a  TEA  CO,  laser  amplifier  requires  a  peak  power 

o 

of  intensity  greater  than  100  kW/cm".  In  our  new  system  we  use  a  3-electrode  I  LA  C  O., 
discharge  tube  which  has  a  higher  gain  than  the  pin-resistor  TEA  CO.,  discharge 
tube.  With  this  improved  system  we  have  obtained  mode-locked  pulses  of  peak  [lowers 
in  excess  of  300  kW . 

(ii)  A  germanium  acousto-optic  modulator  at  Brewster  angle  is  now  being  used.  This 
eliminates  the  power-density  limitation  that  arose  in  the  old  system,  which  used  an 
antireflection- coated  germanium  modulator. 

(iii)  An  antireflection- coated  GaAs  crystal  is  now  placed  outside,  rather  than  inside, 
the  optical  cavity,  which  makes  alignment  of  the  cavity  less  critical.  Also,  since  the 
antireflection  coating  of  the  GaAs  crystal  was  slightly  damaged,  there  is  less  distor¬ 
tion  of  the  TEM  mode  when  the  crystal  is  outside  the  cavity. 

(iv)  In  the  old  system  the  high-voltage  supply  to  the  cavity -dumping  GaAs  crystal 
was  falsely  triggered  by  electrical  noise  generated  by  the  laser  discharge  tube.  In  the 
new  system  a  laser- induced  spark  gap1  is  used  to  eliminate  this  problem. 

(v)  Previously,  continuous  RF  powder  was  supplied  to  the  germanium  mode-locking 
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crystal.  Tins  caused  the  Germanium  crystal  to  heat  up,  thereby  changing  its  acoustic 
resonance  frequencies.  We  are  now  using  pulsed  RF  power  to  minimize  this  problem. 

Lxpe rimental  Arrangement 

Figure  Vl-22  is  a  photograph  of  the  new  ns  TLA  CO,,  laser  pulse-producing  system, 
and  Fig.  Vl-23  is  a  schematic  diagram  of  the  experimental  arrangement.  With  a  tie  volt¬ 
age  of  5  k  V  across  the  C.aAs  electro-optic  modulator  a  fraction  of  the  energy  from  the 
beam  can  be  switched  out  at  the  Brewster-angle  germanium  plate  (CieP,  see  Fig.  Vl-23). 
This  fraction  of  the  beam  was  guided  by  two  flat  mirrors  to  a  gold-coated  minor  with 
2-m  ratlins  of  curvature,  which  focuses  the  beam  inside  the  3-electrode  laser  amplifier, 
flic  tube  is  operated  between  10  Torr  anti  400  Torn.  A  fraction  of  the  beam,  before 
going  into  the  amplifier,  is  reflected  out  from  a  beam  splitter  and  detected  by  a  copper- 
doped  liquid -helium -cooled  detector  with  a  rise  time  of  <1  11s.  The  fraction  of  the  beam 
that  has  passed  through  the  amplifier  is  detected  by  a  gold -doped  liquid-nitrogen-cooled 
detector  with,  a  rise  tune  of  ~1  11s. 

Experiments 

a.  Generation  of  Nanosecond  TLA  CO.,  Laser1  l’ufses 

l  be  laser  cavity  has  an  optical  length  of  approximately  1.  88  m.  To  achieve  forced 
mode  locking.  RF  driving  power  at  2.  5  ms  with  4  W  peak  power  was  supplied  to  the  ger¬ 
manium  acousto-optic  modulator.  The  3-electrode  discharge  tube  inside  the  cavity  is  set 
to  trigger  at  2  ms  after  the  RF  power  is  on.  We  have  obtained  mode-locked  pulses  <2  ns 
wide  (FW11M)  with  a  peak  power  of  ~500  kW.  Figure  VI -24  shows  a  typical  train  of 
mode-locked  pulses  detected  by  a  copper-doped  detector  and  displayed  on  a  Tek¬ 
tronix  7  904  oscilloscope. 


Fig.  Vl-24. 

Typical  mode -locked  pulses  from  the  laser 
oscillator  (5  ns/div). 


The  switching  out  of  individual  pulses  is  accomplished  by  using  the  laser-induced 
spark  gap,  which  is  normally  filled  with  prepurified  nitrogen  at  100  psi.  The  coaxial 
cable  was  charged  up  to  15  kV.  We  can  vary  the  temporal  triggering  of  tins  gap  In 
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altering  the  nitrogen  pressure,  the  distance  between  tin*  two  electrode  of  tin  s|mi  u  .  up, 
or  but1..  Will’ll  tin-  spark  nap  is  triggered,  a  -q  i.iiv  voltage  puls.  ol  7.  5  kV,  twice  tin 
length  of  cable,  is  produced,  l -  .  •  is  supplied  to  the  GaAs  crystal,  which  rotates 
tin-  polarization  ol'  tin-  desired  number  ol'  mode-locked  pulses.  These  mode -lo.  .;>■  I 
pulses  are  then  reflected  out  ol  tile  train  at  lilt  Mi .  w  stcr-angle  german.  nn  plate. 

b.  Xmphfication  ol  Nanosecond  I  f.  \  1'00  Laser  Pulses 

Three  pulses  from  a  train  of  mode-locked  pulses  were  -witched  out  and  curled 
through  an  amplifier  tube.  These  pulses  me  measured  before  and  after  they  propagated 
through  the  laser  amplifier.  figure  Vl-25  shows  input  and  output  signal.-  detected  b\ 


copper-doped  and  gold-doped  germanium  detectors,  r.spetti  ■  lv.  1  hi  output  signal 
from  the  gold-doped  germanium  detector  w  is  rtifieially  delayed  200  ns  by  us.  ol'  a 
coaxial  delay  line.  By  using  the  add-mode  on  .  lektronix  454A  oscilloscope,  w.  w.  r. 
hie  to  display  both  signals  on  .  single  trace. 

In  our  preliminary  studies  we  have  found  that  th.  first  puls,  of  our  three  1...  .  i  pul-, 
i  of  sufficient  power  density  to  saturate  the  I  .sei  implifier.  I  urther  work  on  this 
aspect  of  th.  experiment  will  be  deserin.  .i  in  ,  future  report. 
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V.  Manichaikul 
Lxpe r intent 

We  have  previously  reported  on  the  generation  and  amplification  of  high-intensity 
nanosecond  pulses.'  Two  or  three  of  these  pulses  were  produced.  I  hey  were  trom 
the  1  >(  1  h )  transition,  2  ns  wide  (FW1IM),  separated  by  12  ns.  When  these  pulses  were 
focused  into  a  three-electrode  laser  amplifier  as  shown  in  I*  if*.  \1-12,  a  peak  intensity 
of  2-3  MW/cm‘  was  obtained.  A  beam  splitter  was  used  so  that  the  intensity  of  the 
pulses  could  be  monitored.  The  input  and  output  detectors  were  as  shown  in  Kin.  \  1-12. 
In  this  experiment  the  detected  input  signals  were  delayed  100  ns  by  using  00  It  ol 
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Pig.  V 1- 1 2.  Lxperimcntal  arrangement  for  amplification  of 
high-intensity  ns  pulses.  (See  V.  Manichaikul.') 


lxG-8  cable.  The  add  mode  of  a  Tektronix  oscilloscope  was  used  to  display  the  signals 
for  both  input  and  output  pulses  on  the  same  screen.  Hie  two  detectors  were  calibrated 
against  each  other  by  comparing  the  oscilloscope  picture  of  the  input  and  output  pulses 
without  discharge  exciting  the  three-electrode  laser  amplifier.  Figure  \1-I3a  shows 
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I  ig.  \  1-13.  Output  and  input  of  AHN  experiment.  Total  pressure: 

200  Terr.  Gas  mixture:  CO,  :  N, :  He  =  X :  4 : 1 0  0 . 

Intensity  of  input  pulse:  0.  75  (MW/em2)/div.  Inten¬ 
sity  of  output  pulse:  1.12  (MW/cm  ")/div.  Time: 
50  ns/div.  (a)  Amplifier  off.  (h)  Amplifier  on. 


the  oscilloscope  display  for  this  ease. 

In  order  to  probe  the  temporal  evolution  of  the  three -electrode  laser  amplifier,  we 
first  fired  the  amplifier  and,  after  a  chosen  delay  time,  the  oscillator.  In  general,  the 
oscillator  was  fired  ~30  ps  after  the  onset  of  the  discharge  in  the  amplifier  for  two  rea¬ 
sons:  first,  we  wished  to  avoid  the  effects  on  our  measurements  of  the  shock  waves 
generated  by  the  discharge.  Second,  we  wished  to  be  certain  that  the  symmetric 
stretching  (SS)  and  bending  (B)  modes  of  CO,  had  equilibrated  with  each  other  at 
slightly  above  the  kinetic  temperature  of  the  gas.2 

Measurements  on  the  amplification  of  high -intensity  ns  pulses  were  made  at  200  dorr 
of  C02:N2:He  mixtures.  The  ratio  of  these  mixtures  was  CO,:N,:He  =  X:4:  100,  where 
t  02  partial  pressure  was  varied  from  3.  5  to  35  Torr  partial  pressure.  Small-signal 
gain  of  this  three -electrode  laser  amplifier  in  each  case  was  measured  by  a  cw  CO, 
laser.  2 

Results 

figure  A  1-1 3b  illustrates  the  input  and  output  pulses  when  the  amplifier  is  turned 
on.  hour  such  measurements  were  made  and  their  average  was  taken  at  each  CO 
partial  pressure  studied.  We  have  found  that  the  RG-8  cable  used  for  the  time  delay 
introduces  some  distortion  in  the  input  signals.  This  distortion  can  be  accounted  for 
it  the  first  (second)  pulses  of  the  input  and  output  pulses  from  the  amplification  mea¬ 
surements  are  compared  with  the  first  (second)  pulses  of  the  input  and  output  pulses 
when  the  amplifier  was  evacuated. 
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1  i^.  Vl-14.  I  .SG.,/ LSGj  vs  partial  CO, 
pressure. 


Table  \rI-l.  Experimental  results. 


No. 

r 

pco2 

(Torr) 

SSG 

±10% 

LSGj 

±10% 

Ta 

vv 

(“lx) 

AN 

u 

N  (0) 
u 

1 

3.  5 

0.  50 

0.  30 

1180 

0.  1  2 

2 

6.  5 

0.  70 

0.  36 

1150 

0.  11 

3 

11.0 

1.  50 

0.  76 

1155 

0.  14 

4 

20.  0 

1.00* 

0.  52 

845 

0.  11 

5 

23.  0 

1.  50* 

0.75 

860 

0.  13 

6 

24.  0 

2.  00* 

0.96 

890 

0.  14 

7 

27. 0 

2.  50* 

1.15 

935 

0.  15 

S 

35.  0 

2.  00 

0.  81 

830 

0.  10 

Notes:  'Not  measured  directly:  calculated  from  LSGj  and  the 
peak  intensity  of  the  pulse. 


SSG  =  ^-j — —  is  the  small-signal  gain  across  the  tube, 
in  n 

Here  the  intensity  is  less  than  1  W/cm". 

LSGj,  large-signal  gain  of  the  first  pulse. 

Ta  ,  temperature  of  the  asymmetric  stretching  mode  cal- 
vv  culated  from  SSG. 

AN 

- — ,  fractional  depletion  of  the  00"!  population  by  an 

N^iO)  ns  pulse. 

AN  ,  calculated  from  the  large-signal  gain  and  the  inten- 
U  sity  of  the  pulse. 

N:u(0),  obtained  from  SSG. 
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l  igure  VI-14  shows  LSG^LSGj  vs  the  partial  pressure  of  CO  studied.  We  ha 


ve 


LSG.  =  i  ~  i  _ 
in  i 


where  i  1,2,  with  1  and  2  representing  first  and  second  pulses.  LSG  (LSG  )  is 
the  large-signal  gain  of  the  first  (second)  pulses.  The  following  observations  can  be 
made  from  these  measurements,  (i)  LSG.,/ LSG j  is  less  than  unity.  This  is  to  be 
expected,  since  the  first  pulse  had  depleted  a  fraction  of  the  population  from  the  00°  I 
level  of  COz.  (n)  The  ratio  LSG^LSGj  is  approximately  0.  8  at  COz  partial  pressure 
of  3.  5  Torr  and  the  ratio  increases  slowly  to  0.  9  as  CO.,  partial  pressure  increases  ' 
to  20  Torr  or  higher,  which  is  as  expected,  since  the  00°  1  level  of  CO  was  being 

re  populated  by  the  higher  O  n  m  levels  at  a  rate3  that  is  directly  proportional  to  the 
C02  partial  pressure. 

Table  VI-1  gives  other  experimental  results  of  interest.  We  found  that  the  large- 
signal  gain  of  our  pulses  is  approximately  one -half  the  small -signal  gain,  and  the  frac¬ 
tional  depletion  of  the  00°  1  level,  - 2L,is  between  0.  10  and  0  15 

N  (0)  ‘  ’ 

A  theoretical  model  for  the  amplification  of  high -intensity  nanosecond  pulses  is  being 

developed.  We  shall  present  the  theory,  and  make  a  comparison  of  theory  and  experiment 
in  a  future  report. 
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Recent  advances  in  high-pressure  gas  discharge  technology  have  made  possible  the 
deposition  of  as  much  as  300  joules/liter  in  carbon  dioxide  laser  media.  Because  of 
the  several  vibration-rotation  degrees  of  freedom  of  the  CO,  molecule,  this  energy  is 
stored  in  a  great  many  molecular  states,  and  therefore  efficient  extraction  ol  it  requires 
optical  pulse  lengths  that  are  large  compared  with  the  kinetic  collision  times  governing 
the  energy-exchange  processes  in  the  medium.  Operation  of  high-pressure  devices  as 
oscillators  yields  as  much  as  BO  joules/liter  from  the  afterglow  of  a  pulsed  discharge, 
and  quasi  cw  operation  has  given  hundreds  oT  joules/liter  in  1 0 -1 00  ps  pulses.  Extrac¬ 
tion  in  ns  pulses  is  far  less  efficient.  In  this  report  we  summarize  theoretical  results 
from  a  multitemperature  kinetic  model  formulated  to  describe  ns  pulse  amplification  by 
devices  operating  at  pressures  above  1  atm.  Numerical  results  for  1 -atm  and  5-ntm 
pulse  amplifiers  are  presented. 

Amplification  of  pulses  comparable  to,  or  faster  than,  kinetic  collision  times 
requires  consideration  of  the  polarization  of  the  molecules,  and  of  inertial  effects  in  the? 
molecular  dipoles  (for  example,  see  llopf  and  Rhodes  ).  1  heoretical  models  incorpo¬ 

rating  only  two  vibration  states  and  the  full  rotation  spectrum  will  be  adequate 
descriptions.  If  energy  is  to  be  extracted  efficiently,  the  pulse  length  must  be  several 
collision  times,  in  which  ease  the  coherence  effects  can  be  neglected,  and  the  medium 
can  be  described  by  a  rate -equation  model. 

In  the  rate-equation  limit  the  growth  of  <  plane  wave  in  a  transversely  uniform 
medium  with  non  resonant  loss  o  can  be  described  by  a  first  -order  differential  equation 
in  distance 

=  -d  +  ffM  Nc  jloo0!,,!’  |  -  Y  [10°0.  J  ij-  I.  (1) 


C.  W.  Werner 
U.  Wildman 
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where  Hit*  stimulated  emission  cross  section  is 


yp  '  ( U— (x  )  +  v*- 

o  1  o 

with  Nc  the  C02  density,  Rj  the  degeneracy  of  the*  Jth  rotational  state 

^  =  2J+1,  (3) 

and  (n]n2n3,,j|  is  the  fractional  CO,  population  in  the  state  of  these  quantum  numbers 
(that  is,  the  diagonal  element  of  the  density  matrix  for  a  single  molecule).  The  spon¬ 
taneous  emission  time  written  here  is  that  for  the  entire  band;  an  individual  line  has 
a  matrix  element  proportional  to  the  rotational  matrix  element 

•>  J  for  a  P(.l)  line  (J'  =  ,1  -  ]) 

T  j  j ,  =  < 

J  +  l  for  an  li(.J)  line  J'=,|  +  n 

U'e  shall  neglect  frequency  pulling  effects,  although  in  high-gain  systems  they  will  be 
important  if  the  input  pulse  is  detuned  appreciably. 

Ue  assu-i  e  that  the  molecular  kinetics  can  be  described  adequately  by  rate  equations 
in  which  only  binary  collisions  are  important.  Kven  with  the  rate-equation  model,  the 
six  degrees  of  freedom  of  the  CO,  molecule  in  its  electronic  ground  state  give  rise  to 
so  many  important  vibration-rotation  states  that  the  problem  would  be  intractable  with¬ 
out  further  simplifying  assumptions.  The  fact  that  the  molecule  is  reasonably  harmonic 
in  its  low  vibrational  states,  and  that  the.  interaction  of  vibration  and  rotation  is  weak 
allows  us  to  treat  the  relaxation  of  the.  various  degrees  of  freedom  substantially 
independently.  U'e.  also  make  use  of  the'  observation  of  Osipov  and  Stupochenko2  that 
relaxation  of  molecular  vibrations  from  a  nonequilibrium  distribution  takes  place  in  two 
phases:  first,  a  rapid  relaxation  to  quasi  equilibrium  in  which  the  various  normal 
modes  of  the  mole-cule  acquire  a  Bolt/mann  distribution  of  excitation,  which  can  corre¬ 
spond  to  a  temperature-  far  different  from  the  kinetic  temperature  of  the  gas  and  second, 

a  slow  relaxation  of  these  quasi -equilibrium  distributions  to  the-  kinetic  tempera¬ 
ture. 

Since  we  are  concerned  with  amplifiers  in  which  the  pumping  takes  place  over  a  time 
scab-  that  is  large  compared  with  the  kinetic  collision  times,  it  is  reasonable  to  assume 
that  prior  to  the  arrival  of  the-  electromagnetic  pulse  the  vibrational  states  are  dis¬ 
tributed  according  to  the  partial  equilibrium  distribution 


i  i  --  n, 

lnj  n2n3  J  =  (l-s)(l-b)  ( 1  —a )  s  1  (n.,  +  ] ) 


QPR  No.  Ill 


100 


//  ^  - 


(VI.  APPLIED  PLASMA  RESEARCH) 

where  s,  b,  and  a  are  Boltzmann  factors  for  the  symmetric  stretch,  bending,  and 
asymmetric  stretch  modes  ( v  ^ ,  v9,  v^)  of  the  CO,  molecule.  We  have  taken  the  state 
[n,n.,n3J  to  include  all  of  the  states  [n^n^n^],  of  which  there  are  n9  +  1  (f  represents 
an  angular  momentum  around  the  symmetry  axis  of  the  molecule,  and  hence  can  take 
on  the  values  -n9,  -n„  +  2,  -n  +4,  ....  n9-2,  n  ). 


s  =  exp  "  lTT 

B  s 


b  =  exp 


kB1b 


a  =  exp  ~  E7T 

B  a 


The  rotation  states  also  reflect  a  Boltzmann  distribution 

e  ,1  f  £  /  hcB  \  r  hcB  T/  T , ,  ,1 

_V2n3’  Jj  =  n2n3 Jz  8j  J(J+1  )J  • 


where  we  assume  that  the  rotational  constant  B  is  independent  of  the  vibrational  state 
(C09  has  cB  =  11606  MHz  and  1  1698  MHz  in  the  upper  and  lower  laser  levels,  respec¬ 
tively).  Doppler  broadening  of  the  laser  lines  is  less  than  1%  of  the  homogeneous  line- 
width  at  1  atm,  so  the  velocity  distribution  of  the  molecules  will  be  neglected. 

Passage  of  an  optical  pulse  will  introduce  deviations  from  these  distributions.  In 
particular,  a  fast  pulse  will  create  a  "hole"  in  the  state  [00°1,J'J  (that  is,  it  will 
depress  the  population  below  that  given  by  (5)),  and  a  "peak"  in  1 1 0°  0 ,  Jj  because  of  the 
stimulated  emission  process.  Judicious  approximations  allow  a  description  of  the  kin¬ 
etics  in  terms  of  variables  giving  the  average  occupations  of  the  three  vibrational 
modes  and  the  depths  of  the  "holes"  in  both  vibration  and  rotation.  In  particular,  we 
assume  that  the  two  laser  states  have  the  forms 


[00°  1 J  =  z”a  +  a 


[lCPOj  =  Z  s  +  p  (11 

and  that  the  other  states  retain  their  previous  occupation  probability  exclusive  of  nor¬ 
malization 


|n]n9n3j  =  7.^  s  (n.ll)  b  "a  \  ( n  3  n9n3)  *  laser  state. 
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The  normalization  condition  requires 
Z'1  -  (l-a-(?)(l-s)(l-b)2  ( 1  -a). 


Defining  the  occupation  fractions  for  the  individual  modes 

x  =  Z  [ninqj  (l4) 

mn 

yn  =  T  [mnq  ]  ( J  5) 

mq 

zni  =  r  lmnciJ-  (16) 

nq 

and  with  tile  assumptions  (10),  (11),  and  (12),  we  find  the  following  expressions. 

Xq  =  U-0-P)U-«)aq,  q  #  0  or  1 

x0  =  1  ~  a  -  o  +  (c+p)a  =  (l-a-p)(l-i)  +  p  ►  (J7) 

x.  =  ( 1  -cz-p)(  1  -a)  a  +  a 
1  J 

yr  -•  (l-G-p)(n+l  )1  r  (18) 


zm  =  (1  -a-p)(l  -s)sm,  m  4  0  or  1 

70  =  1  ~  s  -  p  +  (a+p)s  =  (l-o-p)(l-s)  +  a  >  (]9' 

7-j  =  (l-c-p)(l-s).s  +  p. 

These  distributions  are  illustrated  in  Fig.  VI-8. 

While  this  assumption  ot  "holes"  in  single  vibrational  states  is  a  convenient  approx 
inflation,  the  corresponding  ansatz  lor  the  rotational  distribution  is  supported  experi¬ 
mentally  by  the  work  of  Cheo  and  Abrams.3  They  have  found  that  the  rotational 
relaxation  may  be  J -independent  and  all  rotational  levels  are  thermalized  in  one 
collision  time,  so  that  the  expression 


with 


|00° 1, J'j  =  [00° 1 


«• 

“  7~v  2  R,p  cxp 

+  e  • 

- 

L  \  B  >7  J 

> 

(20) 
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zr‘  =  <*-«>  (i^r)  OI) 

correctly  parametrizes  solutions  of  the  model.  £  measures  the  depth  of  the  "hole" 
in  the  rotational  sublevel  depleted  by  the  radiation;  the  other  levels  are  populated 
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Pig.  VI-8.  Assumed  distributions  of  CO,,  normal  vibrational 
mode  excitations. 


in  proportion  to  a  Boltzmann  distribution  scaled  in  amplitude  by  1 so  that  the  net 
vibratio  lal-state  population  is  held  constant  as  £  varies.  A  similar  expression  is 
assumed  for  the  lower  level: 


[  1  0°  0,  J  J  =  [io°o 


2  ffj  oxp 


he  B  \ 

kBTrj 


J(J+1) 


(22) 


with 


Z 


-1 

r 


(i  -n) 


hcB  \ 

k„T  ' 
B  r  / 


This  rotational  distribution  is  illustrated  in  Pig.  VI-9. 


(23) 
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J  0  2  4  6 


t'lg.  VI -9.  Assumed  distributions  of  CO  rotational -state 
excitations.  ^  1 

The  Landau- Teller  assumption4  ,ha,  the  dependence  ol  energy  exchange  cross  sec 
fens  18  that  of  harmonic  oscillator  matrix  elements  can  be  used  to  determine  the  cros- 

unTtlcV  1  Pr°CCSSOS  “°  °rdUrS  ,rom  th»  We  have  taken  i„t« 

account  the  following  processes. 

^  _rl'  in  the  i>_  modn : 


(n]n2n3)  +  M^=±(n].n;,-l,n  )+M  +  e. 


Intermode  V-V  between  r,  and  vv: 


an 

(n]t  n2+I,  n3)  +  N2(v)  ^ - (n^n^  +  N2(v+1)  +  18  cm_1 


Intermode  V-V  between  v3  and  3i^: 


ab 


(nln2'  n3  +  1)  +  M  ^ — —  (nrn2+3,n3)  +  M  +  346  cm-1 
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Intramode  V-\  in  i ■  : 


k 

an 

< n j n 2 ,  n3  +  l)  +  (n^n!,,  n^)  ^^(n^r^)  +  (njni,,  n^  +  1 ) 


R-T: 


k 

r 

(j)  +  M  (j')  +  M. 

In  addition,  the  lower  laser  state  has  been  assumed  to  have  a  V-V  relaxation  of  inde¬ 
terminate  nature  which  lias  been  modeled  by  a  simple  exponential  decay.  Other  V-T 
processes  could  be  included,  but  the  principal  loss  rate  from  both  r,  and  n  is  bv  the 

J  i\  J 

V'2  mode.  Also,  because  of  the  close  coupling  of  the  v  and  modes  by  the  I’crmi  reso- 
nance  (that  is,  large  cross  sections  for  the  conversion  of  one  member  of  a  Fermi  reso¬ 
nant  pair  into  the  other),  we  have  assumed  that  the  i ^  and  »■,  vibrational  temperatures 
are  equal.  The  proper  variable  for  the  description  of  the  combined  bath  of  states  is 


Q  =  2  S  +  B, 


(24) 


where  S  and  B  are  the  average  occupancies  of  r  and  : 


s  =  r 

m 

mz 

m 

-  7  -1 

s 

V 

1  - 

b  =  r 

n 

ny  = 
J  n 

y _1  _ 

2b 

v  1 

-  b  ’ 

A  =  r 

n 

qx  = 

M  q 

7"1  - 

a 

v  1 

-  0 

+  a . 


(25) 

(26) 

(27) 


For  simplicity,  we  have  also  assumed  =  2  so  that 


s  =  bc. 


(28) 


)  he  derivation  ot  the  equations  for  A,  N,  Q,  a,  f3,  £,  q  and  the  kinetic- rotational  energy 
Per  particle  is  straightforward  but  tedious.  Neglecting  the  a  +  (3  terms  in  ( 1 7 ) -( 1 9 )  com¬ 


pared  to  1,  we  find  (oj=uo): 


M  =  _K 

at'  an 


A(  N+l )  -  ( A+l  )N  t 


-e  /k  T 
an  B 


-R 


ab 


1  +  imr]  'A 


(i+ ')  -'"1 >(?)" 


3  -t..,./k„T 

e 


-R. 


29) 
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r  -  +3ltm>[1  +  j]T3j  *(§♦  i)3  -  c~C“b/kB' 


-my  - 


dt7  =  Hna  A(N+I)-(A+1)N 


'£an/kBT 


8t'  ~  “2Kau°  Hi 


3P 

at'  ~  _2UbbP  +  Hi 


=  -YrS  -  (1-1)  IL/lOG0 1 


at'  “  "Yrn  +  d-n)  R./[i<r  oj 

(I  +  +  ^n)  a?  (kBT)  =  t  e  ii  /jLlJi 


b> 


+  4jceanHan  A(N+1 )  ~  (A+l  )N  e 


+  ij<  e  .  R  ,  1  + 

c  ab  ab  \  B  +  2 


A(f+i)  "  (A+J)(f)3  0 


where 


V4nn  =  VL,  " 


R  ,  =  I  <k.k^M^P 

n/r  M  ab 


=  i  *MkLM,p 


H  =  k  P 

a  a  aa 


B  =  B(Tb=T) 
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Ki  =  tr<uo) 


(42) 


Equations  1  and  29-36  (which  are  in  the  canonical  form  of  a  set  of  hyperbolic  equa¬ 
tions  if  t1  is  the  retarded  time  t  -  z/c)  have  been  solved  numerically  for  a  30%  CO, 
gas  mixture  at  1  atm  and  5  atm  total  pressure.  The  input  pulse  was  a  1  ns  (FWHM), 

1  MW/cm  Gaussian  shape.  Initial  conditions  were  calculated  under  the  assumption 

of  equilibration  of  Q  at  T  ,  and  of  T  =  T  sufficient  to  give  the  small -signal  gains 

o  ci  n 


Fig.  VI-10.  Peak  intensity  of  amplified  short  pulse  in  CO, 

media  of  known  small -signal  gain  at  1  and  5  atm 
total  pressure.  Gas  mixture:  0.30:0.05:0.69 
(CO,:N,:He).  Input  pulse:  Gaussian  shape,  1  ns 

C  C.  / 

FWHM,  peak  intensity  10  W/cm  . 
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illustrated;  these  are  typical  of  those  to  be  expected  in  various  high-pressure 
discharges.  Peak  pulse  intensity  as  a  function  of  depth  in  the  amplifier  is 
shown  in  Eig.  VI-10,  and  the  output  pulse  shapes  are  shown  in  Pig.  VI-11. 


Fig.  VI-11.  Temporal  pulse  profiles  for  amplifiers  of 
Fig.  VI-10,  4%/cm  small -signal  gain  at 
1  atm,  3% /cm  at  5  atm. 


Note  that  there  is  substantial  broadening  of  the  pulse  at  1  atm,  while  at  5  atm  the  output 
is  a  reasonably  faithful  diplicate  of  the  input,  even  after  amplification  by  more  than  1000 
in  intensity.  Note  also  the  extremely  large  saturation  intensity  at  5  atm:  elementary 
considerations  of  the  saturation  process  show  that  it  should  scale  approximately  as 
the  square  of  the  total  pressure. 
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H.  A.  Haus,  "A  Theory  of  Forced  Mode  Locking" 


A  THEORY  OF  FORCED  MODE  LOCKING* 
by  H.  A.  Haus** 

Abstract: 

A  theory  of  forced  mode  locking  is  set  up  on  the  basis  of 
an  expansion  of  the  electromagnetic  fields  in  the  optical  cavity 
in  terms  of  the  cavity  modes.  Mode  locking  is  viewed  as  injection 
locking  of  the  cavity  modes  via  generation  of  frequency  sidebands 
by  the  mode  locking  element.  The  approach  leads  naturally  to 
simplifying  assumptions  that  are  legitimate  in  most  practical 
cases.  In  the  case  of  a  homogeneously  broadened  laser  medium,  a 
simple  differential  equation  is  obtained  for  the  mode  locked  pulse 
train  for  an  arbitrary  periodic  mode  locking  modulation.  The 
analogy  of  the  equation  with  well, known  solved  problems  of  quan¬ 
tum  mechanics  is  used  to  determine  the  nature  of  the  solutions. 

For  any  given  modulation,  usually  more  than  one  mode  locked  solu¬ 
tion  is  found.  In  the  case  of  sinusoidal  modulation,  one  of  the 
solutions  reduces  to  that  previously  derived  by  Siegmann  and 
Kuizenga.  The  higher  order  solutions  are  shown  to  be  unstable. 

The  case  of  mode  locking  of  an  inhomogeneous ly  broadened  laser 
medium  is  solved. 

*  Work  supported  by  Joint  Services  Electronics  Program  (Contract 
DAAB07-71-C-0300) .  US  Army  Research  Office  -  Durham  (Contract 
DAHCO  4-72-C-00  44 )  . 

**  Electrical  Engineering  Department  and  Research  Laboratory  of 
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Introduction 


The  theoretical  investigations  of  mode  locking [see  Ref.  1 
for  a  large  list  of  references]  may  be  roughly  divided  into  two 
groups,  (a)  those  which  treat  the  mode  locking  process  in  the 
time  domain,  e.g.  the  paper  by  Kuizenga  and  Siegman[2],  and 
(b)  those  which  treat  the  phenomenon  in  the  frequency  domain 
(notably  the  paper  by  McDuff  and  Harris [3]).  In  this  paper  we 
present  a  useful  extension  of  the  mode  locking  analysis  in  the 
frequency  domain.  When  setting  up  the  analysis  of  saturable 
absorber  mode  locking  by  a  somewhat  novel  approach  we  soon  dis¬ 
covered  that  this  approach  lent  itself  very  well  to  a  better  ’ 
understanding  of  the  forced  mode  locking  phenomenon  and  permitted 
the  solution  of  previously  unsolved  problems  of  mode  locking 
theory. 

This  paper  is  devoted  to  forged  mode  locking.  In  particular 
the  following  issues  are  addressed:  1.  Mode  locking  of  a  homo¬ 
geneously  broadened  laser  by  sinusoidal  amplitude  modulation. 

Well  known  results  are  obtained,  but  this  part  of  the  paper  is 

intended  mainly  to  show  the  correspondence  with  previous  work. 

* 

2.  Mode  locking  by  sinusoidal  amplitude  modulation  of  an 

inhomogeneously  broadened  laser.  3.  Modulation,  other  than 

sinusoidal,  of  a  homogeneously  broadened  laser.  4.  Stability 

of  the  mode  locked  pulses  of  a  homogeneously  broadened  laser. 

« 

The  presentation  of  a  theory  of  mode  locking  from  a  dif¬ 
ferent  point  of  view  needs  some  justification  in  view  of  the 


large  literature  on  the  subject.  It  is  believed  that  the  approach 
presented  here  offers  the  following  advantages: 

(a)  The  theory  of  coupled  oscillators  is  a  well  known  one 
and  offers  insights  which  are  not  as  easily  obtained  using  the 
approach  of  pulse  propagation  through  nonlinear  media. 

(b)  The  theory  of  injection  locking  has  been  recently  dis¬ 
cussed  very  lucidly[4].  The  analysis  of  the  stability  of  locked 
oscillators  can  be  transferred  almost  unchanged  to  the  treatment 
of  the  stability  of  mode  locked  pulse  trains. 

(c)  Simplifying  assumptions  concerning  the  character  of  the 
laser  medium  are  suggested  by  the  new  approach  which  lead  to  very 
simple  differential  equations.  The  very  simplicity  of  these  dif¬ 
ferential  equations  then  allows  the  treatment  of  problems  not 
hitherto  attempted.  The  present  paper  illustrates  two  such  ap¬ 
plications:  mode  locking  of  an  inhomogeneous ly  broadened  laser 
and  stability  analysis  of  the  mode  locked  pulse  train. 
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I.  Mode  Locking  as  Injection  Locking 

Injection  locking  of  a  cavity  with  a  single  mode  containing 
a  negative  conductance  may  be  analyzed  with  the  aid  of  the  equiv¬ 
alent  circuit  of  Fig.  1.  Here  the  admittance  y^(|v|)  represents 
the  active  admittance  of  the  gain  element  which  is  in  general 
complex  and  a  function  of  the  magnitude  of  the  complex  voltage 
V,  and  Yc  is  the  cavity  admittance,  a  function  of  u.  If  a 
current  source  Ig  of  frequency  w  is  applied  across  the  circuit 
as  indicated,  then  the  equation  for  this  injection  locked  oscillator 
is  given  by 

t 

Is  -  lYc(w)  +  Yl(u,  | V| ) ]  V  (1.1) 

Equation  (1.1)  can  be  adapted  to  the  case  of  a  mode  locked 
optical  cavity  involving  the  interaction  among  many  axial  modes. 

The  identification  of  equivalent  circuit  elements  of  each  reso- 
ance  circuit  (representing  one  axial  laser  mode)  with  laser  and  • 

optical  cavity  parameters  is  carried  out  in  Appendix  I. 

The  electric  field  in  the  cavity  is  expanded  in  the  set  of 
axial  modes  of  the  empty  cavity,  with  mode  pattern  e^(r);  the 
expansion  coefficients  become  the  voltages  of  the  different 

resonance  circuits.  One  assumes  at  the  outset  that  mode  locking 
is  successful  and  that  all  modes  oscillate  near  their  respective 
resonance  frequencies.  The  frequency  separation  of  the  oscilla¬ 
tions  is  dictated  by  the  injection  signals  which  are  by  definition 
separated  by  wM,  the  frequency  of  the  mode  locking  element. 

l 

■ea*  -- 


* 


Whether  or  not  the  solutions  thus  found  are  physically  acceptable, 
i.e.  stable,  and  the  starting  assumption  of  successful  mode 
locking  is  indeed  justified,  must  be  left  for  a  stability  analysis 
which  is  carried  out  at  the  end  of  the  paper. 

These  assumptions  imply  that  the  electric  field  in  the  cavity 
can  be  expanded  as  follows: 

-  -  jkwMt 

E(r,  t)  *  e  °  E  Vk  e  "  ek(r)  +  c.c.  (1.2) 

k 

where  the  k-th  mode  has  the  frequency 

wk  =  ^  WM  +  wo*  (1.3) 


jw  t 

The  fast  time  dependence  exp  °  has  been  explicity  extracted 
where  u>0  is  to  be  identified  with  the  frequency  of  the  central 
mode  (k  =  0) ,  the  mode  nearest  to  the  peak  of  the  laser  gain 
profile  (peak  of  the  negative  conductance  Re  Y_ )  . 

Xj 

Each  complex  amplitude  obeys  a  set  of  equations  that 

can  be  summarized  in  a  parallel  equivalent  circuit.  The  equiva¬ 
lent  circuits  describing  different  modes  are  coupled  by  the  side¬ 
band  generation  of  the  mode  locking  element,  i.e.  the  driving 
current  for  the  k-th  cavity  mode,  I£k'  contains  contributions 
of  V^(i  ?  k) .  We  shall  disregard  couplings  as  caused  by  other 
effects.  Thus,  in  particular,  we  shall  assume  that  the  gain  of 
the  laser  medium  (represented  by  Y  )  is  time  independent  and 
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therefore  no  sidebands  are  produced  on  frequency  components 

amplified  by  the  laser  medium.  This  in  turn  implies  that  the 

«* 

population  relaxation  time  is  long  compared  with  the  cavity 

transit  time  so  that  the  laser  medium,  while  affected. by  the 
time  average  optical  power  in  the  cavity,  does  not  respond  to 
instantaneous  power  fluctuations.  (This  assumption  was  also  made 
in  Ref.  2.)  Generalizing  (1.1)  to  the  k-th  cavity  mode  with 
anjiitude  one  has: 

|yCk<“k>  +  W1  Vk  =  rsk  U-4> 

Consider  first  the  dependence  of  Y„  (w,  )  upon  io.  .  The 

CiC  K  k 

k-th  cavity  mode  has  a  quality  factor  Q  (assumed  to  be  the 
same  for  all  modes)  and  a  resonance  frequency 

wko  *  k  Au  +  u0  *  (1.5) 

where  Aio  is  the  separation  between  the  resonances  of  the 
(empty-axial)  cavity  modes.  Since  the  frequency  of  excitation 
is  given  by  (1.3),  one  has 

t 

-  Au 

— -  .  (1.6) 

“o 

Next,  consider  of  (1.4).  If  the  cavity  is  mode  locked 

by  the  time  dependent  loss  of  a  mode  locking  element,  then 


Sk 


m(t)  E  V.  e 
l  * 


(1.7) 


k-th  Fourier  component 


where  m(t)  has  the  time  dependence  of  the  loss  modulation  (with 
period  2tt/wm)  and  its  magnitude  is  a  measure  of  the  "loading" 

of  the  cavity  by  the  mode  locking  element.  (See  Appendix  I  for 
details.) 

In  (1.7)  all  modes  A  affect  the  mode  k  equally  (m(t) 
is  k  and  £  independent)  an  assumption  justified  if  the  mode 
locking  element  is  near  one  of  the  mirrors.  Other  positions  of 
the  mode  locking  element  must  be  handled  differently  as  shown  in 
Section  V.  The  expression  (1.7)  may  also  be  written, 


*Sk  " 


Gp  E  M(k  ~ 

C  l  1 


(1.8) 


where  M  ^  is  the  (k  -  i) - th  Fourier  component  in  the 

Fourier  decomposition  of  m(t); 


(k)  7k  wMt 


m(t)  =  E  Mvk;  e~  M 
k 


(1.9) 


Combining  (1.4),  (1.6),  and  (1.7)  we  obtain 


1  +  k  j  2Q 


wM  -  Aw 


w. 


Vk  +  V“k>  Vk 


,  jiuMt 

Gc  [m(  t)  E  Vt  e  M  ] 

I  k-th  component 


(1.10) 


mmurn 


*»  w’UfttiJiiy' 
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Finally,  consider  the  admittance  of  the  laser  medium  (for 
details,  see  Appendix  I)  .  The  laser  gain  is  frequency  dependent, 
and  a  Lorentzian  line  shape  will  be  assumed.  The  admittance  may 
be  normalized  to  G_  so  that 


W 


-  G- 


1  +  j 


“k  "  wo 


Wt 


(1.11) 


where  «L  is  the  laser  medium  linewidth,  and  g  is  the  magnitude 
of  the  negative  conductance  at  line  center  normalized  to  Gc> 

We  have  assumed  that  the  center  frequency  of  the  laser  line 
coincides  with  one  of  the  cavity  mode  oscillations.  If  the  mode 
spectrum  is  dense  as  it  must  be  to  get  good  mode  locking,  this  is 
not  a  serious  restriction . 

We  shall  now  make  the  most  important  assumption  of  this  study, 
which,  in  fact,  is  the  assumption  crucial  for  obtainment  of  the 
simple  equations  of  this  treatise.  We  assume  that  (1.11)  can  be 

expanded  in  wk  -  w0/wL  and  only  terms  up  to  second  order  are 
retained 


W  *  -  Gcs 


(1.12) 


This  assumption  implies  that  the  pulse  widths  predicted  by  the 
ensuing  analysis  must  be  long  compared  with  1/^.  This  assumption 
is  usually  justified  in  practice. 
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The  conductance  g  is  power  dependent.  If  one  normalizes 
the  field  pattern  e^  appropriately,  the  power  P  on  the  travel¬ 
ling  waves  of  the  (standing  wave)  modes  propagating  in  one  axial 

direction  can  be  set  equal  to 


(1.13) 


for  a  homogeneously  broadened  medium;  is  the  saturation  power 

of  the  laser  medium  in  general  a  function  of  if  the  pulse 

spectrum  is  narrow,  as  assumed,  the  w  -dependence  can  be  disregarded 
(1.12)  is  correct  strictly  speaking  only  if  all  modes  act  on  the 
laser  medium  in  an  identical  way,  i.e.  if  the  medium  is  concentrated 
near  one  of  the  mirrors.  The  phenomenon  of  "spatial  hoi  a  burning' 
is  thus  eccluded .  Assumption  (1.12)  may  not  be  too  restrictive 
even  for  a  laser  medium  occupying  a  larger  portion  of  the  cavity 
length . 

Suppose  at  first  that 

m(t)  =  2M [1  -  cos  wM  t]  (1.14) 

the  modulation  is  purely  sinusoidal.  Then 

XSk  =  GC  MlVk  -  1  -  2  Vk  +  Vk  +  1]*  <1-15) 

Adhering  to  this  example,  and  introducing  (1.15),  (1.3),  (1.5)  and 

(1.12)  into  (1.10)  we  obtain  the  second  order  difference  equation: 


II.  Solution  of  Sinusoidal  Mode  Locking  of 
Homogeneously  Broadened  Laser 

We  make  a  further  approximation  so  as  to  enable  us  to  find 
simple  solutions  to  the  difference  equation  (1.16).  If  many 
modes  are  to  be  locked,  and  if  the  variation  of  the  gain  over 
a  frequency  range  corresponding  to  the  mode  spacing  is  small, 
then  the  difference  equation  can  be  replaced  by  a  second  order 
differential  equation.  One  has 


Vk  +  1  "  2Vk  +  vk  -  l 


where  Vk  is  now  replaced  by  a  continuous  function  of  the  variable 

frequency  d0Viation  from  “o'  and  is  the  frequency 

separation  of  the  spectral  lines.  The  modes  are  identified  by 

the  continous  variable  u  =  lim  k„,  where  k  is  ^ 

k-*» 

number,  which  assumes  w  •♦■o 

M 

both  positive  and  negative  values.  One  obtains  for  (1.16) [5]: 


1  -  g 


(i  ^  («  +  g) 

'  WL  /  WL 


-  (6  +  g)  v  =  Mw  2 
L  dw 


(2.1) 


where  6  is  a  measure  of  cavity  "detuning",  proportional  to  the 

difference  between  the  modulation  frequency  and  cavity  resonance 
frequency  spacing: 


Eq.  (2.1)  is  particularly  simple,  if  we  require  that 


6  +  g  =  0. 


This  implies  that 


1  Vo 

-  g  - 

20  .  UL 


(2.3) 


(2.4) 


The  modulation  frequency  deviates  from  the  cavity  frequency  by 
an  amount  proportional  to  the  laser  medium  dispersion.  In  fact, 
(2.4)  means  that  the  modulation  frequency  is  adjusted  to  equal  the 
cavity  mode  frequency  separation  as  modified  by  the  laser  medium. 
Hence  (2.4)  may  be  viewed  as  excitation  by  a  modulation  tuned 
to  the  actual  mode  separation,  i.e.  resonant  modulation.  With 
(2.3)  satisfied,  (2.1)  has  a  pure  real  coefficient  on  the  right 
hand  side.  The  resulting  equation  is  identical  with  the  harmonic 
oscillator  equation  of  quantum  mechanics,  where  V  plays  the 
role  of  the  wave  function  and  w  the  role  of  the  spacial  vari¬ 
able. 


d2  V 
du2 


(2.5) 


We  can  now  bring  to  bear  the  entire  formalism  of  the  quantum 
mechanical  harmonic  oscillator  on  the  problem  of  loss  modulation 
in  a  homogeneously  broadened  medium.  The  solutions  of  (2.5)  are 


13 


the  well  known  Hermite  gaussians[6] 


V(w)  =  Hn  \  exp  - 

W/ 

where 


The  role  of  the  energy  eigenvalues  En  of  the  harmonic  oscillator 
is  played  by  the  parameter  g  -  1  in  (2.5).  We  find  for  the 
quantization  of  this  parameter 

2  /Fig  wM  (n  +  1/2) 

E  «  g  -  1  -  - - -  (2.8) 

“L 

Equation  (2.6)  defines  "supermodes"  or  "mode  locking  modes" 
of  the  mode  locked  oscillator.  The  higher  the  order  n  of  the 
"supermode",  the  more  structure  in  the  frequency  spectrum.  Since 
the  Fourier  transform  of  the  spectrum  (2.6)  leads  to  the  same 
time  dependence,  we  find  that  mode  locked  solutions  may  exist  of 
higher  and  higher  temporal  structure  exhibiting  more  and  more 
pulses.  The  lowest  order  supermode  was  obtained  by  Siegman  and 
Kuizenga(2]  and  the  higher  order  supermodes  by  Haken  and 
Pauthier[7]  and  others[8,9J.  We  sf all  show  in  Section  VI  that 
all  higher  order  supermodes  (n  >  0)  cannot  be  excited  in  the 


steady  state  because  they  are  unstable.  Equation  (2.8)  determines 
the  power  in  the  mode  locked  cavity.  Indeed  by  solving  for  g 
from  (2.8)  and  then,  using  the  dependence  of  g  upon  power  P, 
(1.13) ,  we  find  the  following  relationship  for  P: 


P_  1  +  y  +  /2y  +  y2 
Li 


where  y  =  —  (2n  +  1)  2  M 
2 


In  the  limit  of  a  small  modulation  coefficient  M,  (2.9)  becomes 


w. 


(n  +  1/2) 


(2.10) 


Let  us  now  return  to  an  interpretation  of  the  present 
solution  in  terms  of  the  physical  model  used  as  a  starting  point 
of  the  investigation ,  the  injection  locking  model.  We  must  re¬ 
member  that  the  differential  equation  was  used  as  a  mathematical 
approximation  to  the  difference  equation,  and  that  the  solution 
obtained  must  be  interpreted  as  an  approximate  solution  of  the 
difference  equation  if  the  replacement  is  made:  V(to)  -*■  and 

the  solution  (2.6)  is  assigned  values  only  at  discrete  frequencies 
(o^,.  Mode  locking  has  led  to  simultaneous  oscillation  of  many  modes 
in  a  homogeneously  broadened  system,  even  though  the  unmodulated 
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system  would  have  oscillated  only  in  a  single  mode.  The  simul¬ 
taneous  oscillation  in  many  modes  is  allowed  because  of  injection 
locking.  The  injection  currents  are 

*Sk  ■  GCM  \  -  1  -  2  Vk  +  Vk  ♦  1  • 

» 

Figure  2  shows  a  plot  of  the  real  part  of  the  two  sides  of  (1.16) 
vs  w.  (the  imaginary  part  vanishes  because  of  assumed  resonant 

Iv 

modulation),  showing  explicitly  the  discrepancies  RelY^w^)  + 

Y  (w.  ) ]  which  are  balanced  by  the  modulation-produced  "source" 

L  k 

current  normalized  to  the  oscillation  voltage,  wv 

The  results  obtained  thus  far  are  well  known.  One  can  extend 
our  expressions  to  cover  reactive  (phase)  rather  than  loss  (ampli¬ 
tude)  modulation  by  replacing  M  by  an  imaginary  quantity ,  again 
arriving  at  results  derived  previously.  The  mathematical  steps 
necessary  for  such  modifications  are  interesting,  but  are  not 
presented  here  because  no  new  results  are  obtained [10] .  Instead, 
we  shall  concentrate  here  on  problems  whose  solutions  were  not 


obtained  previously. 
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III.  Inhomogeneous ly  Broadened  Medium 


If  the  medium  is  inhomogeneous ly  broadened,  then  the  power 
in  any  particular  cavity  mode  determines  the  pulldown  of  the  gain 
line  at  that  particular  cavity  frequency  (provided,  of  course, 
that  the  separation  of  the  cavity  modes  is  large  compared  with 
the  homogeneous  linewidth  of  the  medium) .  The  adaptation  of 

(2.1)  to  the  case  of  an  inhomogeneous ly  broadened  medium  is  ac¬ 
complished  by  changing  the  character  of  the  dependence  of  the 
negative  conductance  g  upon  intensity.  To  first  order,  the 
reactive  contribution  is  unaffected  by  saturation.  The  equation 


(3.1) 


where  PL  is  a  conveniently  defined  saturation  intensity.  Eq. 

(3.1)  is  a  nonlinear  differential  equation  for  the  voltage  V. 

It  is  a  difficult  equation  to  solve;  yet  it  is  possible  to  gain 
considerable  insight  into  the  nature  of  its  solutions  by  giving 
it  a  physical  interpretation.  Let  us  separate  the  complex  voltage 
amplitude  into  an  amplitude  and  a  phase  factor 


V  =  A  e 


j* 


(3.2) 


■  * 


Taking  the  second  derivative  of  (3.2)  one  obtains 


=  A  ej*  +  2j i  k  *  j*  A  -  i*  A  el* 

dw2 


(3.3) 


Introducing  (3.2)  and  (3.3)  into  (3.1)  one  has 

..2  1  v 

gD  1  -  — 

..  .  wM2 

M  «M2  (A  -  <£2  A)  =  1 - p— . - . -  A 

ll  +  — 


(3.4) 


M  uM2  (aJ  +  2A  i)  »  Ji  U  +  g  )  A 


(3.5) 


Here  we  have  written  (  )  for  d/dw  in  order  to  emphasize 
the  analogy  with  the  equations  of  motion  of  a  particle  in  polar 
coordinates;  the  frequency  variable  plays  the  role  of  "time",  the 
amplitude  A  is  analogous  to  "radius".  The  force  field  is 
"radius"  and  "time"  dependent.  If  one  assumes  that  the  driving 
frequency  is  matched  to  the  frequency  separation  of  the  cavity 
modes  as  modified  by  the  medium,  then  6  +  gQ  =  0  and  the  force 
field  becomes  a  central  one.  In  this  case,  one  may  define  a 
potential  function,  the  derivative  of  which  gives  the  force.  The 
potential  function  plotted  against  "radius"  A  as  a  function  of 
"time"  03  is  shown  in  Figure  3.  The  potential  hill  travels  as 
a  function  of  time  towards  the  origin  and  then  the  origin  becomes 


a  hill  in  its  own  right.  If  one  looks  for  a  mode  locked  "super- 
Iriode" ,  one  looks  for  the  motion  of  a  particle  that  starts  with 
zero  velocity  at  a  time  t  =  0  (zero  slope  at  center  frequency) , 
and  then  approaches  the  origin  (A  =  0)  where  it  comes  to  a  full 
stop. 

Consider  first  the  case  in  the  absence  of  modulation.  In 
the  particle  motion  analog,  this  corresponds  to  a  massless  parti¬ 
cle.  The  particle  "looks"  for  a  position  in  which  it  is  exposed 
to  zero  force  field.  The  particle  stays  at  the  top  of  the  hill 
starting  infinitesimally  to  the  left  and  then  searches  for  the 
position  in  which  the  force  is  infinitesimally  small  yet  directed 
from  right  to  left.  When  the  potential  hill  has  made  it  to  the 
origin,  the  particle  has  followed  it  there  and  then  stays  there 
forever.  This  case  can  be  solved  very  simply  analytically  by 
setting  the  righthand  side  of  equation  (3.4)  =  0. 

Consider  next  the  case  when  there  is  modulation,  at  first 
the  case  of  very  small  mass  (modulation)  .  Near  t  =  0,  when 
the  force  field  is  strong,  the  particle  seeks  out  positions  in 
the  potential  field  at  which  it  is  exposed  to  a  small  force  to 
the  left.  This  means  that  the  particle  starts  to  the  left  of  the 
hill  at  t  =  0,  i.e.  at  a  smaller  radius  than  the  massless  par¬ 
ticle.  In  the  language  of  mode  locking,  this  corresponds  to  a 
power  of  the  mode  near  line  center  that  is  smaller  than  the  power 
in  the  absence  of  mode  locking.  Near  the  origin,  the  particle 
has  acquired  a  finite  velocity  and  kinetic  energy.  In  order  to 
come  to  a  full  stop  at  the  origin,  it  has  to  expend  the  kinetic 


■  . . . 
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energy  by  climbing  up  a  hill.  This  means  that,  before  the  particle 
comes  to  full  rest,  the  origin  must  have  already  risen  above  its 
surroundings.  Translated  into  the  language  of  the  mode  locked 
"mode"  this  means  that  the  frequency  bandwidth  of  the  mode  is 
wider  than  in  the  absence  of  modulation  (in  the  language  of  the 
equivalent  particle  the  origin  of  the  radial  coordinate  develops 
into  a  hill  before  the  massless  particle  has  reached  the  origin) . 

The  model  is  useful  also  for  determining  what  happens  when 
the  strength  of  mode  locking  is  increased  excessively.  This 
corresponds  to  a  very  massive  particle.  The  potential,  as  defined, 
does  not  involve  the  mode  locking  amplitude  M  and  hence,  in  the 
equivalent  language  of  the  particle,  the  initial  hill  becomes 
smaller,  the  heavier  the  particle.  A  supermode,  symmetric  at 
line  center,  corresponds  to  an  initial  condition  for  the  analog 
starting  at  rest  to  the  left  of  the  potential  hill.  A 
particle  which  starts  initially  at  rest  slightly  to  the  left  of 
the  potential  hill  may  be  too  massive  to  make  it  near  the  origin 
before  the  potential  distribution  has  changed  beyond  recall, 
preventing  the  particle  from  acquiring  sufficient  kinetic  energy. 
The  particle  never  makes  it  to  the  origin.  No  solution  exists 
for  excessively  high  values  of  M. 

There  are  also  higher  order  solutions,  if  the  modulation  is 
not  too  strong.  One  may  start  the  particle  on  the  down  slope  of 
the  hill.  (Note  the  peak  intensity  at  center  frequency  for  this 
mode  is  less  than  for  the  lowest  order  supermode.)  The  particle 
may  go  through  the  origin  and  climb  up  the  hill  on  the  other  side. 
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After  oscillating  back  and  forth,  the  particle  can  be  brought  to 
rest  at  the  origin.  Such  solutions  exist  only  for  particles  of 
finite  mass  (i.e.  finite  modulation) . 

Equation  (3.4)  has  two  adjustable  parameters,  gQ  -  1, 
the  excess  small  signal  gain  at  line  center  (w  =  0)  and  the 
normalized  modulation  parameter 


M  WM 


9o  UL 


For  a  range  of  these 
two  parameters,  solutions  for  the  particle  motion  were  obtained. 
Figure  4  shows  plots  of  V»^  vs  normalized  frequency  w/wM,  for 
various  choices  of  the  parameters  /  — —  ~  and  gn  -  1.  Figure 


g  u. 

^o  L 


5  shows  the  range  in  the  (gQ  -  1) , 
~M  WM 


~  plane  over  which  mode  locked  solutions  are  obtained. 

:  %  wl 

Figure  6  shows  the  decrease  of  the  oscillation  amplitude  at  line 
center  A(0)//P 7,  as  a  function  of  normalized  modulation,  for 

Xj 

different  excess  small  signal  gains. 


IV.  General  Time  Dependent  Loss  Modulation 


The  generalization  to  a  general,  time  dependent,  loss 
modulation  calls  for  re-examination  of  the  "injection  current" 
expression  (1.7) .  This  expression  suggests  that  an  analysis  of  the 
mode  locking  problem  in  the  time  domain  may  lead  to  simplifications. 
Indeed,  if  one  defines 

iS(t>  5  1  *sk  e  «•!> 

one  finds 

Ig(t)  =  -  Gc  m(t)  v ( t)  .  (4.2) 

3  iw^t 

where  v(t)  »  E  vc  e  is  the  time  dependent  amplitude  of 

i 

the  electric  field  near  one  of  the  mirrors.  Equation  (4.2)  states 
simply  that  the  time  dependent  source  current  is  the  product  of 
the  modulated  time  dependent  conductance  Gc  m(t)  and  the  field 
amplitude. 

All  of  this  is  straight  forward  and  would  not  profit  us 
very  much  if  one  could  not  transform  the  left  hand  side  of  (1.4) 
into  the  time  domain  and  obtain  a  simple  result.  In  general,  this 
is  not  possible.  However,  with  the  approximations  already  made 
f°r  an^  YL^wk^  we  ma^  carry  out  the  transformation 

and  obtain  a  very  simple  result.  Using  the  same  approximations 
that  led  to  (1.6),  one  obtains  in  the  time  domain: 
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1  -  g 

i  +  -L- 

+  (g  +  5) 

_d_ 

v  ( t)  =  m(t)v(t) 

(4.3) 

II  “l2  dt2 

dt 

In  the  case  of  "resonant”  excitation,  g  +  5  =  0,  one  obtains 
a  differential  equation  in  time  for  v(t)  with  a  simple  quantum 
mechanical  analog,  A  particle  with  the  mass  2m/iT2  =  w  2/g 
moves  in  the  potential  well  m{t)  where  the  role  of  the  spatial 
variable  x  of  quantum  mechanics  is  now  played  by  the  time.  The 
sinusoidal  modulation  problem  of  the  homogeneously  broadened  laser 
could  have  been  solved  in  the  time  domain,  with  1  -  cos  oo.,t-+  4(oj,,t)  2 
and  the  hermite  gaussian  solutions  would  have  been  obtained, 
another  way  of  observing  the  invariance  of  hermite  gaussians  with 
regard  to  Fourier  transformation.  Note,  however,  that  a  solution 
of  the  inhomogeneous ly  broadened  laser  mode  locking  problem  would 
have  been  difficult  in  the  time  domain. 

Solutions  (4.3)  can  be  obtained  once  one  assumes  special 
forms  of  the  modulation  function.  Consider  for  example  a  very 
deep  modulation  of  the  mode  locking  crystal  which  shuts  off  trans¬ 
mission  entirely  outside  an  interval  of  length  t.  This  is  shown 
in  Figure  7.  Equation  (4.3)  for  the  square  potential  well  is  the 
quantum  mechanical  problem  of  an  electron  in  a  perfectly  reflecting 
box.  The  eigenfunctions  and  eigenvalues  of  this  problem  are  well 
known  and  are 


V(t)  =  cos 


mr 


t 


(4.4) 


(WWPW'J'  '•'  •i?r wr  »9*  m  n,wup  nw**rV*'! r—mmu^  ,  /  vrr^ * — *1^1. ^  a.j  i„  ■  i * 1 1  jju  lu'HfB.iu*  ■ 


in  the  time  domain.  The  eigenvalues  are 
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(4.5) 


In  the  time  domain  the  mode  locked  pulse  is  entirely  confined  to 
the  time  period  of  length  x. 

If  the  square  wave  modulation  is  modified  to  one  of  lesser 
depth,  one  has  the  problem  of  an  electron  in  an  open  square  well 
with  well  known  eigenfunctions (6 ] . 

One  may  think  of  many  other  modulation  forms  for  all  of  which 
it  is  easy  to  arrive  at  qualitative  conclusions  as  to  the  shape 
of  the  mode  locking  pulse  and  the  possibility  of  mode  locking. 

The  viewing  of  the  mode  locking  problem  in  the  time  domain 
affords  other  insights.  Consider  for  example  the  question  as  to 
what  approximations  were  made  when  the  discrete  spectrum  of  the 
cavity  oscillations  was  replaced  by  a  continuum  and  the  modulation 
of  loss,  which  was  sinusoidal,  was  replaced  by  the  operator  d2 /dw2 . 
In^the  time  domain,  this  operator  is  replaced  by  the  function 
~  2  *  •  If  one  considers  the  periodic  problem,  one  has  to  make 
two  modifications:  (a)  one  must  return  to  the  full  time  dependence 
.of  the  modulation  and  replace,  in  the  time  domain,  the  parabola  by 
the  cosine  function.  (b)  One  must  look  for  solutions  of  the  time 
dependent  problem  which  are  periodic. 

Clearly  the  introduction  of  a  periodic  modulation  function 
in  lieu  of  the  parabolic  one,  introduces  a  periodic  well  in  time 
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space.  The  problem  is  now  a  quantum  mechanical  problem  in  a 
periodic  potential.  If  the  modulation  is  strong  enough  so  that 
the  eigenvalues  lie  near  the  bottom  of  the  well,  one  obtains  a 
set  of  solutions  in  each  of  the  well  bottoms,  which  is  almost  the 
same  gaussian  solution  that  has  been  obtained  previously.  The 
eigenfunction  consists  of  isolated  pulses  in  time,  repeating 
periodically.  On  the  other  hand  when  the  modulation  decreases, 
the  walls  of  the  periodic  wells  decrease  and  for  a  given  strength 
of  modulation,  solutions  are  found  with  overlap  between  the  wave 
functions  within  each  of  the  periodic  wells.  The  mode  locking 
is  not  strong  enough  to  suppress  the  amplitude  between  pulses. 
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V.  General  Position  of  Mode  Locking  Element 

When  the  mode  locking  element  is  moved  a  distance  zq  away 
from  one  mirror,  where  Wp  zq/c  is  comparable  to  one  or  larger, 
one  may  not  assume  that  the  modes  of  interest  have  roughly  the 
same  field  patterns  thoughout  the  mode  locking  element. 

If  the  time  varying  conductivity  of  the  mode  locking  element 
is  crM(r,  t)  ,  the  injected  current  Igk  is  proportional  to 

_  jfcWw1 

-  I  /  aM(r,  t)  e^*  •  e^  dv  e 

k-th  Fourier  component 

and  m(t)  of  (1.7)  has  to  be  replaced  by  an  expression  proportional 
to 

i\A(t)  «  /  oM(r,  t)  e^  •  ek*  dv 

which  depends  on  k  and  l.  Now 
(k  Aw  +  oj  ) 

e.  «  sin  - 2—  z 

c 

and  using  a  corresponding  expression  for  e^,  one  finds  that 
ek*  *  ei  contains  some  rapidly  spatially  varying  terms,  which 
integrate  out  to  zero  and  a  more  slowly  varying  term 
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Spatial  integration  over  the  "short"  mode  locking  element  replaces 
z  by  zq ,  where  zQ  is  the  position  of  the  element  measured 
from  one  mirror.  The  injection  current  ig(t)  in  the  time  domain, 
as  defined  by  (1.8),  becomes 


The  time  dependent  injection  signals  are  produced  by  the  delayed 

and  advanced  versions  of  the  mode  locking  modulation,  a  very 

obvious  result.  What  is  remarkable,  however,  is  that  we  may  now 

effect  a  simple  generalization  of  the  results  of  the  preceding 

section  to  obtain  the  mode  locking  solutions  for  this  problem. 

Suppose  we  have  a  square  wave  modulation  as  shown  in  Figure  8.  If 

A  Zn 

the  time  delay  (advance)  —  —  is  one  quarter  of  the  period 

WM  c 

of  uM  (center  of  the  cavity,  Aw  =  wM)  the  number  of  potential 
wells  has  doubled  as  shown  in  Figure  8  but  their  depth  has  de¬ 
creased.  The  number  of  mode  locked  pulses  doubles  provided  the 
depth  of  the  wells  is  still  sufficient  to  allow  "trapped"  solutions. 
If  the  modulation  is  of  a  more  general  shape,  displacement  of  the 
element  causes  variations  in  the  shape  of  the  effective  potential 
well,  leading  to  distortions  in  the  wave  function  (pulse  shape) . 

The  present  approach  can  be  used  to  ascertain  pulse  shapes  in  a 
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VI.  Stability  of  Solutions  of  Homogeneously  Broadened  Laser  Medium 

In  this  section  we  study  the  stability  of  the  supermodes 
with  respect  to  arbitrary  perturbations.  We  shall  follow  the 
analysis  of  stability  of  injection  locking  of  Kurokawa[4] .  In 
this  analysis,  the  assumed  steady  state  is  perturbed,  the  per¬ 
turbation  is  taken  to  be  a  slow  function  of  time  (compared  with 
the  time  variation  of  the  unperturbed  oscillation) . 

An  instantaneous  frequency  and  growth  or  decay  rate  is  ob¬ 
tained  for  the  perturbation  and  the  transient  time  evolution  of 
the  perturbation  is  studied  for  net  growth  or  decay.  If  it  can 
be  shown  that  there  is  no  net  growth  the  injection  locking  is 
stable,  otherwise  it  is  unstable. 

We  shall  extend  this  analysis  to  the  study  of  the  stability 
of  "steady  state"  mode  locking  solutions.  Here  all  oscillations 
are  coupled  to  each  other  through  the  injection  locking  process. 

The  stability  study  is  concerned  with  the  time  evolution  of 
these  coupled  perturbations. 

Study  of  transients  implies  the  replacement  of  the  frequencies 
^k  of  inoc*es  by  +  where  J?k  is  complex;  -Im(n^)  >  0 
implies  growth  of  excitation.  Further  the  "instantaneous"  per¬ 
turbation  frequency  itself  is  a  function  of  time.  Thus  even 

temporary  growth  may  be  offset  by  eventual  decay. 

The  objective  then  is  to  obtain  an  equation  for  a  complex 
perturbation  <5Vk  of  the  steady  state  Vk  and  the  frequency 
nk  of  the  perturbation  which  is,  in  general,  a  function  of  the 


phase  4>  between  6v  and  V..  The  perturbation  obeys  the 
same  difference  equation  (1.16)  as  except  that  ^k  must 

be  replaced  everywhere  by  tok  +  flk  and  g  may  be  perturbed 
because  an  amplitude  perturbation  changes  P  and  thus  changes 
g.  All  a>k  dependent  parameters  may  be  expanded  to  first  order 
in  the  «k's  if  the  flk's  are  assumed  small  (small  growth 
rates  and  frequency  deviations) . 

We  shall  disregard  all  derivatives  of  the  parameters  of  (1.4) 
with  respect  to  ftk,  except  that  of  Im  .  This  is  tan¬ 

tamount  to  disregarding  all  energy  storages  other  than  the  elec¬ 
tromagnetic  energy  storage  of  the  cavity  modes.  If  the  oscillation 
of  the  perturbation  6Vk  of  the  k-th  axial  mode  occurs  with 
the  frequency  wk  +  then  the  admittance  to  be  associated  with 

«Vk  is 
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*ck  <“k  +  V  -  yck<“k>  +  i  nk  - 

“k 

or  in  the  continuum  limit 

Y_(to  +  (1)  =  Yr  (oj)  +  jfJ  —  v(6.1) 

wo 

If  we  assume  that  the  unperturbed  solution  is  the  n-th 
supermode,  V^fw)  •  excited  at  synchronism  (i.e.  g  +  6  =  0) 
expand  the  difference  equation  as  indicated  and  go  to  the  continuum 
limit,  we  obtain  the  equation: 


I  M  ^  l  vj  J' 


+  (6g)  /l - -  -  j  JS.  \  V  =  j«  6V 

\  “L  “l  /  wo 


(6.2) 


Next,  consider  the  perturbation  vn<59* 


V« 


V„  « 


i  +  -L 


V  g  / 

n^o  .  /  P 

i  ♦  -lV  k 


(6.3) 


If  the  voltage  Vn(w)  is  so  normalized  that  /|Vn(w)j2  dw  is 
equal  to  the  power,  then 


6P  -  /|Vn(u)  +  6V(w)  |  2  dw  -  /|Vn(w)|2  du  = 


and  thus  we  obtain  for  (6.2) 


2  /  Vn(u))  Re  6V  du 


2  d2  _  ,  gw2  ] ...  /  _ 


...  \  Vo  /  V  (u)Re  6Vd  u 
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It  is  convenient  to  test  first  the  stability  of  a  general 
n-th  supermode  with  respect  to  a  perturbation  6V(w)  proportional 
to  the  m-th  supermode,  6V(U)  =  am  um(oj)  ,  and  study  the  fre¬ 
quency  fi(w)  of  the  perturbation;  here  u^tw)  is  a  properly 
normalized  eigenfunction  of  the  defining  equation  of  the  supermodes 


M  w  2  — -  g  — u  =  -  E  u 

M  dw2  w  2  m  m  r 

Li 


(6.6) 


and  the  eigenvalues  Em  are  given  by  (2.8).  Since  the  eigen¬ 


functions  are  orthogonal,  we  may  set 


/  U  udu>  =  <5 
_  n  m  nm 


(6.7) 


Introducing  6V  =  am  um(uj)  in  (6.5)  and  using  the  orthogonality 
condition,  we  obtain  for  the  frequency  deviation: 


fi  *=  j  —  (E  -  E  ) 

2q  m  n 


(6.8) 


For  ^m  <  En'  i»e*  i°r  m  <  n,  one  obtains  Imft  <  0  and  one 
finds  growth  of  the  perturbation.  Thus  any  n-th  supermode  is 
unstable  with  respect  to  a  perturbation  having  a  dependence  upon 
cavity  mode  number  (w)  corresponding  to  a  lower  order  supermode. 
Hence  all  modes  n  >  0  are  unstable. 

The  only  possible  stable  mode  is  the  0-th  order  (gaussian) 
mode.  In  order  to  prove  its  stability  we  must  show  that  a  general 


r 


J.' |||  s'- jgup 
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perturbation  5V  =  I  a^  u^ (td)  does  not  lead  to  exponential 


m 


growth.  We  find  from  (6.5),  using  V  (w)  =  /p“  u  (w) 

o  o 


'  S  (E* '  E°’  a^M  '  777  2Rea°  f1  -  77  -  J  V“> 


(d 


2  —  jfl  Z  a  u  (a)) 
m  m  m 


(6.9) 


Through  multiplication  by  um((d)  ,  m  =  0,  1,...,  and  subsequent 
integration  one  obtains  a  set  of  coupled  equations  for  the  coef¬ 


ficients  am«  For  m  =  0  we  obtain 


°9  PL  °°  -2j* 

ao  (1  +  e  °>  ■  jO  at 


0  11  +  — 
P. 


P_ 

> 

L. 


(6.10) 


where 


Coo  5  f  uo2 <“>  <  1  - 


WT 


u‘ »  d«d  =  1  - 1  !t. 

2  Wl’ 


(6.11) 


and  $Q  is  the  phase  angle  of  aQ,  i.e.  the  phase  difference 
between  the  0-th  order  component  of  the  initial  perturbation 
ao  Uo*u)  and  the  steady  state.  Note  that  Cqq  >  0  even  though 
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(6.11)  seems  to  imply  the  possibility  of  negative  values.  We 
recall,  however,  that  we  have  expanded  a  Lorentzian,  an  approxi¬ 
mation  valid  only  when  the  pulse  bandwidth  wp  is  small  compared 

with  toT .  If  the  expansion  had  not  been  made,  (6.11)  would  be 
Li 

positive  definite. 

Equation  (6.10)  is  interpreted  by  a  stability  circle  similar 
to  the  one  discussed  by  Kurokawa  when  treating  stability  of 
injection  locking  of  a  single  oscillator,  Figure  9.  For  an  arbi¬ 
trary  initial  phase  $Q,  the  frequency  deviation  and  rate  of 
decay  are  as  shown  in  the  figure.  For  no  initial  phase  is  growth 
observed.  Suppose  the  initial  phase  of  the  perturbation  $o  >  0 
as  shown  in  Figure  9.  Then  Re  ftQ  >  0  and  the  frequencies  of 
the  perturbations  are  greater  than  the  frequencies  of  the  steady 
state  oscillation.  4>  grows  at  the  same  rate  for  all  modes  and 
hence  the  circle  is  described  with  increasing  time  as  shown.  4> 
varies  as  the  circle  is  traversed.  The  amplitude  of  6 V ( to)  decays 
at  the  rate  Im  ft  which  is  always  positive.  The  perturbation 
has  relaxed  when  the  origin  is  reached.  If  the  initial  phase 
4>o  =  ir/2,  Im(ft)  =  0  and  no  growth  or  decay  occurs,  i.e.  there 
is  no  restoring  force  for  a  perturbation  of  this  kind.  A  quadra¬ 
ture  perturbation  corresponds  to  a  phase  perturbation  of  the 
carrier  frequency.  Hence,  the  mode  locked  pulses  have  no  carrier 
phase  stabilization,  just  like  a  free  running  van  der  Pol  oscil¬ 
lator  [11]  . 

Next  consider  the  equations  obtained  for  the  perturbation 
amplitudes  a^^  and  .  The  equations  for  m  >  2  are  all  given 
by  (6.8)  with  n  =  0  because  the  integral 
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Cmo  5  1  um(w)  uo(w)  fl  “  -^7  “  j  “  ]  d“ 

\  WL  WL  / 

vanishes  for  all  m  ^  0,  1,  2.  Hence  decay  is  predicted  for  all 
perturbations  with  m  >  2.  We  have  for  m  =  1  and  2: 


Cmo  Re  ao  "  ^  am  (6.12) 


This  equation  shows  that  Re  aQ  acts  as  a  source  of  these  per¬ 
turbations.  In  fact ,  a  proper  interpretation  of  (6,12)  is  obtained 
only  after  one  recalls  that  jfl  is  equivalent  to  a  time  derivative 
in  the  spirit  of  the  present  analysis  (and  Kurokawa's  stability 
analysis)  which  deals  with  the  time  evolution  of  an  instantaneous 
amplitude  and  frequency. 

Since  the  natural  "frequencies"  of  the  system  of  (6.12)  are 
decaying  exponentials,  (6.12)  predicts  solely  that  the  lowest 
order  perturbation  m  =  0  produces  in  =  1  and  m  =  2  pertur— 


W_ 
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VII.  Conclusions 

The  analysis  of  forced  mode  locking  as  a  problem  of  injection 
locking  led  naturally  to  approximations  which  would  have  been  less 
evident  in  a  different  approach.  In  particular,  the  expansion  of 
the  laser  medium  susceptance  in  terms  of  frequency  deviation  intro¬ 
duced,  in  the  time  domain,  first  and  second  derivatives  as  a 
description  of  pulse  distortion  by  the  laser  medium.  This  ap¬ 
proximation  reduced  every  forced  mode  locking  analysis  to  a 
Schroedinger  equation.  Detuning  of  the  mode  locking  drive  intro¬ 
duces  first  derivatives  which  can  be  removed  by  writing  the  func¬ 
tion  v ( t)  as  a  product  of  an  exponential  multiplying  factor  and 
another  time  function .  The  displacement  of  the  mode  locking  crys¬ 
tal  is  represented  simply  as  a  shaping  and  a  repetition  of  the 
"potential  well".  The  inhomogeneously  broadened  mode  locked  laser 
could  be  treated  by  this  new  approach.  The  stability  analysis  of 
the  mode  locked  solution  reduced  to  a  problem  in  essence  treated 
by  Kurokawa. 

It  is  believed  that  the  potential  of  the  present  analysis 
has  been  hardly  tapped  and  that  many  other  issues  of  interest, 
both  in  forced  mode  locking  and  saturable  absorber  mode  locking, 
will  be  analyzable  using  this  new  approach. 
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In  this  Appendix  we  derive  the  circuit  equations  used  in 
the  body  of  the  paper  for  the  mode  locking  analysis.  We  assume 
an  open  cavity  formed  by  curved  mirrors  in  free  space,  portions 
of  which  are  filled  by  the  laser  medium  and  the  mode  locking 
element.  The  electromagnetic  field  in  the  cavity  obeys  Maxwell's 
equations 

V  x  E  =  -  u  —  (1.1) 

°  3 1 


V  X  H  =  En  -  +  J  (1.2) 

°  3t 


In  these  equations  J  represents  the  perturbation  of  the  cavity 
modes  by  thelosses,  tl.a  laser  medium,  and  the  mode  locking  element. 
We  expand  the  electromagnetic  field  in  terms  of  the  cavity  modes 


with  normalized  field  patterns  ek(r) 

and  hj.(r) 

for  the 

elec- 

trie  and  magnetic  field  respectively. 

These  mode 

patterns 

are 

related  by 

V  *  e'k  =  Bk  hCr) 

(1.3) 

V  x  hk  =  8k  h(r) 

(1.4) 

where  $k  is  related  to  the  resonant  frequency  wko 


of  the  cavity. 
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(1.5) 


These  field  patterns  are  divergence  free. 

V  •  ek  =  0  (1.6) 

V  •  hk  =  0  (1.7) 

In  general,  an  expansion  in  modes  having  divergence  is  also 
necessary [12 ] .  Nowever,  the  effect  of  these  modes  is  negligible 
in  optical  cavities.  We  shall  omit  contributions  due  to  such  modes. 
The  field  patterns  are  orthogonal  and  may  be  normalized. 


f  ek  •  e£*  dv  =  /  hk  •  h^*  dv  =  N  6^  (1.8) 

Here  N  is  a  normalization  constant  whose  choice  is  dictated  by 
convenience.  To  be  specific,  one  may  consider  a  Gaussian  mode 
whose  field  pattern  is  proportional  to 


5k(?)  • 


2  /N 
/¥l  w(z) 


sin  B  ^z 


e 


-r  2/  w2 


(1.9) 


where  L  is  the  spacing  between  the  mirrors  and  w  is  the  beam 
diameter  which  is  the  well  known  function  of  position  z  and 
radii  of  curvature  of  the  mirrors[13].  One  may  assume  that  at 
the  mirrors  the  boundary  condition  is  met  that  the  tangential 
electric  field  vanishes.  The  field  is  approximately  transverse 


to  the  direction  of  propagation,  the  polarization  vector  lies 
in  the  surface  of  constant  phase. 

Use  of  the  mode  patterns  effects  a  separation  of  the  time 
dependence  from  the  spatial  dependence.  One  may  write 

-  - 

E (r ,  t)  =  I  vk(t)  ek(r)  e  °  +  C.C.  (1. 10) 

k 

_  _  jeo  t 

H  (r ,  t)  =  E  iu(t)  h.  (r)  e  °  +  C.C.  (1. 11) 

k  k 

where  we  display  only  the  positive  frequency  parts  and  extricate 
explicitly  the  fast  time  dependence  exp  jwQt  corresponding  to 
the  centev  frequency  of  the  laser  line.  We  indicate  all  fre¬ 
quency  components  from  v/hich  the  fast  time  dependence  has  been 
extracted  by  a  bar  under  the  letter.  When  these  expansions  are 
introduced  into  Maxwell's  equations  and  the  ortogonality  condi¬ 
tions  (1.8)  are  utililized,  one  obtains  differential  equations 
in  time  for  the  amplitudes  i^(t)  and  vk(t),  the  equivalent 
currents  and  voltages. 


(1.12) 


e,  *  dv 
k 


(1.13) 


The  equivalent  capacitance  and  inductances  of  (1.12)  and  (1.13) 
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where  wq  is  the  frequency  of  the  mode  nearest  laser-line  center 
and  the  mode  index  k  assumes  positive  and  negative  values. 

If  mode  locking  is  achieved,  the  frequencies  of  the 

Fourier  components  of  the  field  will  be  evenly  spaced  with  the 
frequency  separation  w^,  where  is  the  frequency  of  modula¬ 

tion  of  the  mode  locking  element.  Because  of  the  resonant  nature 
of  (1.12)  and  (1.13),  the  time  dependence  of  the  k-th  cavity 
mode  will  be  essentially  exp  j  w^t  where  is  the  frequency 

of  the  spectral  component  nearest  the  resonance  frequency  w^.0 
of  the  cavity.  This  implies,  on  one  hand,  that  must  be  near 

Aw  and  that 


k  wM  +  uj 
M  o 


(1.16) 


where  toQ  is  the  frequency  of  oscillation  of  the  mode  nearest 


to  laser  line  center.  In  fact,  the  uq's  in  (1.15)  and  (1.16) 


need  not  be  identical;  but  in  the  case  of  cavity  modes  spaced 
closely  compared  to  the  laser  line-width  one  will  always  find 
a  mode  near  line  center  that  oscillates  essentially  at  its  reso¬ 
nance  frequency.  The  equation  of  the  amplitude  of  the  k-th 
cavity  mode,  whose  amplitude  has  the  time  dependence  exp  j  w.  t 
as  explained  above  is; 


j  +  -  /  J  •  e^*  dv 


BkN 


=  0 


(1.17) 


where  the  "susceptance"  of  the  k-th  mode  is 


)  =  (to  +  cu.  )  C,  - 


'C'"V  -  ““o  T  “V  ^  ■  - - --1’  ■  *  2  —  G 

(“o  +  VLk  “o  . 


where 


Q  =  “kock  =  V; 


(1.18) 


with  G^,  to  be  defined  later,  and  wkQ  and  are  approximated 

as  independent  of  mode  number  k.  Here  we  have  denoted  the  Fourier 
component  of  the  last  term  in  (1.17)  at  frequency  to,  by  a  ver- 
tical  bar  followed  by  the  subscript  k.  This  is  the  driving  term 
for  the  k-th  cavity  mode  which  contains  in  addition  to  the  in¬ 
jected  signal  due  to  the  mode  locking  modulation  a  contribution 
of  the  cavity  loss  and  the  laser  medium  gain.  We  assume  that  the 
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j«ofc 

positive  frequency  part  of  J,  with  e  omitted,  J,  may 

be  written  in  the  form 

J(r,  t)  =  [ac ( r)  +  jwQe0  xL(r)  +  0M(r'  J  £(r'  fc)  '  *I,19) 

The  first  term  is  an  equivalent  conductivity  representing  the 
cavity  losses.  The  second  term  includes  the  complex  polarizability 

X  of  the  laser  medium.  The  last  term  is  the  modulated  conduc- 

L 

tivity  produced  by  the  mode  locking  element,  under  the  assumption 
of  resistive  modulation.  A  modulation  of  a  dielectric  constant 
can  be  treated  by  replacement  of  the  conductivity  by  an  imaginary 
quantity.  In  (1.19)  it  has  been  assumed  that  the  cavity  loss  and 
the  laser  medium  susceptibility  are  time  independent  and  hence 
produce  no  modulation  sidebands  of  the  electric  field.  The  instan¬ 
taneous  variation  in  laser  gain  as  produced  by  a  passage  of  a  mode 
locked  pulse  have  been  disregarded.  This  means  that  the  assumption 
has  been  made  that  the  laser  medium  relaxation  time  T^  is  very 
long  and  hence  the  laser  medium  responds  only  to  the  time  indepen¬ 
dent  power  of  the  mode  locked  pulse  train.  The  k— th  Fourier 
component  of  the  driving  term  of  the  k-th  mode  in  (1.16)  assumes 
then  the  following  form 

(t) 

0 

k-th  Foune 
component 


/  J  •  ek(r)  dv 


=  [G„  +  +  Gp  E  mw(t)  V 


(1.20) 
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where  the  various  coefficients  have  been  defined 

G  =  -  /  o  (r)  e  .  e  *  dv  (1.21) 

u  0ON  c 

yl  ~  —  /  XL(r)  e.  •  e  *  dv  (1.22) 

3  N  K  K  . 


GC  mk  (t)  "  f  CTH(?'  t}  V  *  dV 


(1.23) 


If  the  modelocking  element  is  near  one  of  the  mirrors ,  the 
integral  becomes  k-  and  Z-  independent  and  one  may  then  omit 
the  subscripts  k,  Z.  This  has  been  done  in  the  first  part  of  the 
text.  If  the  laser  medium  and  the  equivalent  loss  are  nonuniformly 
distributed  through  the  cavity,  then  coupling  is  produced  between 
the  different  mode  patterns.  This  coupling  leads  to  injection 
signals  caused  by  the  k-th  Fourier  component  into  modes  of  index 
t  ^  k.  Such  injection  signals  excite  the  mode  off  resonance  and 
therefore  produce  a  negligible  effect  and  will  be  disregarded. 

The  conductance  assigned  to  the  mode  locking  element,  on  the  other 
hand,  because  of  its  time  variation,  produces  Fourier  components 
at  frequencies  other  than  the  frequency  of  the  driving  electric 
field,  and  hence  coupling  between  the  modes  due  to  this  conductivity 
°M  n,ust  be  included.  Introducing  (1.18)  through  (1.23)  into  (1.17), 
we  finally  obtain  the  desired  relationship 
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k-th  Fourier  component  (1.24) 


This  is  (1.10)  of  the  text. 

Next  we  study  the  effect  of  gain  saturation  and  the  approxi¬ 
mations  necessary  to  lead  to  manageable  equations.  The  rate 
equation  for  the  population  inversion  for  a  driving  field  near 
the  resonance  frequency  of  the  medium  is 


at  tx 

Here  p  is  the  matrix  element  of  the  laser  levels.  We  assume 
that  all  Fourier  components  of  the  electric  field  lie  sufficiently 
near  the  resonance  of  the  laser  medium  so  that  all  Fourier  com¬ 
ponents  in  E  may  be  assumed  to  have  an  identical  depleting  effect 
on  the  population  inversion.  If  one  assumes  that  the  relaxation 
time  of  the  laser  medium  is  long  enough  so  that  the  laser 

medium  cannot  respond  to  the  time  dependent  components  of  the  mode 
locked  laser  pulses,  one  may  write  for  the  population  inversion 
approximately 


ft 


T2  n 


(1.25) 


- - . . .  ■ 
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n. 


1  -  2< 


nil2 

* 


>Ti  t2 


(1.26) 


where  <>  indicates  a  time  average.  The  time  independent  com¬ 
ponent  of  the  square  of  the  electric  field  may  be  written 


(1.27) 


We  are  interested  in  the  term,  (compare  (1.22)): 


/  XL  ek  •  ek*  dv  «  /  dv  ek 


2T1T2 


(1.28) 


A  simple  evaluation  is  possible  if  the  laser  medium  is  short  and 
positioned  near  one  of  the  mirrors  of  the  optical  cavity.  Using 
the  mode  pattern  (1.9) ,  one  has  for  the  second  term  in  brackets 


of  (1.27) 


/  dv  (£k  •  iR*) (e^  •  e^*) 


2nLwz  l 


(1.29) 
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where  we  have  made  the  assumption  that  the  propagation  constant 
8  is  identical  for  all  modes  of  interest  an  assumption  valid  if 
th®  integration  is  to  be  carried  over  a  volume  short  compared  to 
the  optical  cavity  length  and  situated  near  one  end  of  the  cavity 
the  length  of  the  laser  medium.  With  this  finding  one  may 
write  for  the  integral  (1.27): 


1  XL  5k 


,*  dv  =  /  x£  ek 


•  e, 


dv 


1  - 


s  K|2 


Pr 


(1.30) 


where  the  saturation  power  P  is  defined  by 


3  N 


ir  L  w! 


T1T2 


(1.31) 


When  the  medium  fills  the  entire  cavity,  the  integral  in  (1.28) 

is  somewhat  more  complicated  and  in  fact  what  is  included  by 

carrying  out  the  integral  is  the  effect  of  spatial  hole  burning, 

an  effect  not  to  be  considered  in  the  present  study. 

One  may  generalize  (1.30)  to  be  valid  approximately  also 

for  large  field  intensities  by  replacing  the  factor  1  -  £  V  2 A> 

n  I, 

by  1/(1  +  l  Vn2/Pfa) .  This  is  the  form  used  in  the  bulk  of  the 
paper. 

It  is  convenient  to  normalize  the  "voltage"  components  so 
that  their  mean  square  amplitude  gives  the  power  in  one  of  the 
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Abstract 


A  closed  from  solution  is  obtained  for  node  locking  of  a 


laser  by  a  saturable 
laser  medium  gain  is 
as  a  function  of  time 


absorber  with  instantaneous 
aaaumcd  time  independent. 


response. 
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I.  Introduction 


In  another  publication,  we  have  developed  a  theory  of  forced 
mode  locking  in  the  frequency  domain  (see  appendix  IV) .  With 
some  modifications  of  this  mode  locking  theory,  it  is  possible  to 
obtain  what  we  believe  is  the  first  closed-form  solution  for 

saturable  absorber  mode  locking. 

As  in  the  case  of  forced  mode  locking,  we  treat  mode  locking 
by  a  saturable  absorber  as  a  form  of  injection  locking.  This 
means  that  the  admittance  mismatch  of  the  resonant  circuit  repre¬ 
senting  a  cavity  mode  is  balanced  by  the  injection  locking  current 
produced  through  interaction  of  the  equivalent  cavity  voltage  (i.e. 
electric  field  amplitude)  with  the  saturable  absorber. 

In  equation  form: 

lYCk(“k'  +  V“kn  Vk  *  TSk-  <la) 

Here  Y  is  the  equivalent  admittance  of  the  k-th  axial  cavity 

mode  and  Y  is  the  equivalent  admittance  of  the  laser  medium. 

L 

The  reader  is  referred  to  appendix  IV  for  details  of  the  derivation. 
The  following  approximations  are  made: 

(a)  The  mode  susceptance  is  expressed  as  a  linear  function  of 
frequency  deviation  from  the  mode  resonance  frequency  anc^ 

the  equivalent  conductance  is  assumed  mode  independent 
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YcK> 


=  G, 


1  +  2  j 


on 
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(1.2) 


(b)  The  laser  medium  equivalent  admittance  has  a  Lorentzian  shape 
that  is  approximated  by  an  expansion  in  frequency 


YL*wk*  ~  Gc  g 
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Here,  g  is  power  dependent,  where  P  is  equal  to  the  sum  of  the 
counter  traveling  powers  of  the  standing  waves  of  each  axial  mode. 


g  - 


l  +  — 


(1.4) 


The  saturation  expression  is  approximately  independent  of  frequency 
if  the  frequency  spectrum  of  the  pulse  is  narrow  compared  with  the 
linewidth . 

(c)  All  Fourier  components  of  the  amplitude  of  the  k-th  axial 
mode  other  than  the  one  near  the  cavity  resonance  are  neglected 
(high  Q  approximation) . 

(d)  One  assumes  that  mode  locking  is  successful,  i.e.  that  a  solu¬ 
tion  of  periodic  pulses  does  exist,  so  that  all  Fourier  components 
are  evenly  spaced  at  separations  u>M.  The  frequency  of  the  k-th 
axial  mode  amplitude  is 


wk  = 


k  WM  +  “o 


(1.5) 
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where  u>o  is  the  frequency  of  the  mode  nearest  medium  line  center 
(set  approximately  at  line  center,  an  assumption  legitimate  if 
the  mode  spectrum  is  dense) . 


(e)  The  injection  current  in  (1.1)  is  assumed  to  be  produced 

by  modulation  products  of  the  equivalent  mode  voltages  (electric 
fields) .  The  loss  element  is  assumed  to  interact  equally  with  all 
inodes  (i.e.  the  element  is  thin  and  near  one  of  the  laser  cavity 
mirrors) . 


jiti)Mt 
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(1.6) 


The  separation  of  the  resonance  frequencies  uko  of  the  axial 
modes  is  constant,  Au,  and  thus 


“ko  "  k  Am  +  <J0.  (1.7) 

(1.5)  and  (1.7)  imply  that  the  central  mode,  k  =  0,  oscillates 

a 

at  resonance,  an  assumption  justified  by  the  fact  that  we  obtain 
a  solution. 

Introducing  (1.2),  (1.3),  (1.5),  (1.6),  and  (1.7)  into  (1.1) 
we  obtain 
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k-th  Fourier  component 


The  right  hand  side  of  (1.8)  is  more  conveniently  expressed  in  the 
time  domain,  as  -  m(t)  v(t)  where 


v(t) 


I  V.  e 
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(1.9) 


is  the  time  dependent  amplitude  of  the  standing  wave  electric  field 
near  the  mirror  (next  to  which  the  saturable  absorber  is  situated) . 
One  may  look,  therefore,  for  simplification  of  (1.8)  by  expressing 
the  entire  equation  in  the  time  domain.  Multiplying  the  left  hand 
side  by  exp  j  k  to„t,  adding  over  all  k  and  noting  that  multi- 

Pi 

plication  by  (j  k  u„) n  is  dn/dtn  in  the  time  domain,  one  has 
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dt 
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v  ( t)  =  -  m(t)  v  ( t) 
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where 
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(1.11) 


The  equation  describes  the  change  wrought  on  a  pulse  of  the 
electric  field  by  the  action  of  the  cavity  and  the  laser  medium 
and  sets  this  change  equal  to  the  change  caused  on  the  pulse  by 
the  time  dependent  element  as  expressed  by  m(t)  .  The  term 
1  x  v(t)  is  proportional  to  the  dis torcion-f ree  decrease  of  the 
pulse  caused  by  the  cavity  loss.  The  laser  gain  produces  a  change 
proportional  to 


10) 
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v. 


If  v  is  a  single  pulse  with  two  inflection  points,  the  second 
derivative  is  a  pulse  with  three  extrema.  Remember 
that  d2/dt2  is  the  Fourier  transform  operator  of  -w2 .  Hence 
this  contribution  expresses  the  deficiency  of  gain  at  larger  fre¬ 
quency  deviations  from  line  center  as  compared  to  the  gain  at  the 
center  frequency.  The  pulse  with  the  three  extrema  subtracts  at 
the  center  of  the  original  pulse  (t  -*•  0,  u  -+•  ®)  and  adds  on  the 
wings  ( 1 1 1  •*  u)  -*■  0)  .  The  pulse  passing  through  the  medium 

is  amplified  and  broadened  (in  time) . 

On  the  right  hand  side  appears  the  (negative)  change  produced 
by  the  time  varying  element,  in  our  case  the  saturable  absorber. 
This  element  has  to  compensate  for  the  shape  change  produced  by 
the  laser  medium. 

For  the  purpose  of  relating  to  the  derivation  in  the  next 
section,  note  that  m(t)  as  produced  by  a  saturable  absorber 
may  be  written  as  Q/QA(t)  where  Q  is  the  cavity  quality  factor 
in  the  absence  of  the  absorber  and  l/QA(t)  is  the  time  dependent 
inverse  Q  of  the  saturable  absorber. 
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II.  The  Saturable  Absorber 


The  rate  equation  for  the  population  difference  n  between 
the  lower  and  the  upper  levels  of  the  saturable  absorber  is 


in.  -  _  2  0  Jv(t)Ji 

3t  T.  A  fiw  A 

A  o 


(2.1) 


v(t)  is  so  normalized  that  | v(t) | 2  gives  the  sum  of  the  powers 
in  the  two  counter  traveling  waves  in  the  cavity  added  over  all 
axial  modes.  Here  is  the  relaxation  time  of  the  absorbing 

medium,  n£  is  the  equilibrium  population  difference,  a  is 
the  optical  cross-section  of  the  absorbing  particles  and  A  is 
the  cross-section  of  the  laser  mode.  If  tne  thickness  of  the 
absorbing  medium  is  0^  and  the  length  of  the  cavity  is  L,  one 
may  define  a  Q  associated  with  the  absorbing  medium. 


2  o.  nc  0. 
o  A  A 


where  the  resonance  frequency  of  a  cavity  mode,  is  set  equal 

to  the  laser  line  center  frequency.  This  Q  is  a  function  of 

intensity  and  is  obtained  from  (2.1)  and  (2.2)  as  a  solution  of 

the  differential  equation.  If  the  relaxation  time  T  of  the  ab- 

A 

sorber  is  fast  compared  with  the  rate  of  change  of  the  intensity, 


we  may  assume  that  the  population  difference  is  an  instantaneous 
function  of  intensity,  and  obtain  an  approximate  expression  for 
the  instantaneous  Q  of  the  saturable  absorber 


l.c,A  2°A  ne 

\  |  V  ( t)  1  2  1 

_  1 

\  .  Mali' 

^0 

3 

1 

1  < 

1  O 

PA 

II 

o 

*o 

- 1 

< 

Ok 

1 

H 

_ I 

(2.3) 


where  we  have  defined  the  saturation  power  for  the  absorber  by 
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The  time-independent  part  of  the  Q  of  the  absorber  may  be  in¬ 
corporated  in  the  cavity  loss.  Thus,  for  the  fast  absorber,  we 
obtain  the  differential  equation 
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We  shall  now  look  for  solutions  of  this  equation  corresponding 
to  mode  locked  pulses.  Remember  that  the  frequency  a>M  is  an 
adjustable  parameter,  to  be  determined  from  the  character  of  the 
problem.  It  will  become  clear  that  no  periodic  solutions  exist 
when  g  +  5  ^  0  and  thus  we  may  state  at  the  outset  that  the  pulses 
must  pick  their  repetition  period  so  that  g  +  5  =  0.  The  remaining 
equatior 


can  be  recognized  as  the  equation  of  motion  of  a  particle  of 
displacement  v(t)  in  a  potential  well 


-  -  (1  -  g) 
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v*  +  I  -2_  £ 
4  Qa°  PA 


(2.7) 


If  the  particle  is  launched  at  the  well  height  0,  at 


(2.8) 


with  zero  velocity  it  moves  tc  ci;e  origin  and  stops  there.  This 
solution  is  symmetric  in  time.  A  oulse  like  solution  for  the  motion 
cf  the  particle  is  then 


as  can  be  ascertained  by  substitution  of  (2.9)  in  (2.6). 

This  solution  is  an  isolated  pulse.  A  succession  of  periodic 
pulses,  of  any  desired  period  21/^/271“"-  g)  <  TR  <  -  is  obtaine, 
by  launching  the  particle  at  a  lower  height.  If  xp  <<  TR  = 
then  the  single  pulse  is  an  excellent  approximation  to  one  period 
of  the  periodic  pulse  train.  Well  mode  locked  pulses  with  good 
time  separation  are  of  this  character.  Since 


the  energy  in  the  pulse  is  given  by  2Tp  vq2.  The  power  P  of 

che  laser  is  then  P  =  2 Tp  vo2/TR.  Introducing  this  expression 

into  (2.10) ,  one  obtains  an  equation  between  P  and  r 

lP 


(2.12) 


where  we  have  defined  the  coefficient 


(2.13) 


(2.11)  and  (2  12)  supplemented  by  the  dependence  of  the  negative 
conductance  g  upon  power  P,  (1.4) 


g 


yield  three  equations  for  the  three  unknowns  P,  g,  and  the 
pulse  width  ip. 

From  (2.11),  we  note  that  g  is  less  than  unity.  This  means 
that  the  laser  is  below  threshold  with  respect  to  the  linear  loss 
(loss  :.n  absence  of  laser  power)  .  The  laser  oscillates  because 
the  bleaching  of  the  absorber  reduces  the  loss  below  the  linear 
loss.  The  equivalent  injection  locking  voltages  have  to  be  equal 
to  the  difference  between  the  voltages  produced  by  the  positive 
and  negative  impedances.  Because  this  combination  of  impedances 
is  below  threshold  over  the  entire  mode  spectrum  the  injection 
voltages  are  all  of  the  same  sign.  This  finding  has  to  be  con¬ 
trasted  with  forced  mode  locking,  where  the  center  portion  of  the 
spectrum  is  above  threshold,  the  wings  are  below  threshold  and 
the  equivalent  injection  locking  voltages  change  sign  as  one  pro¬ 
gresses  from  the  center  of  the  line  to  its  edge. 

(2.11)  mav  be  used  to  eliminate  g,  and  from  the  remaining 
two  equations  (2.12)  and  (1.4)  one  may  obtain  two  equations  for 
P/PL  and  o)L  Tp  respectively. 

_L  -L  =  1/(goKl 

gO  PL  (D.  Tp  +  - - - 

"t  TP 

Solutions  of  1/u)tt0  vs  g  <  are  shown  in  Figure  1.  It  is  clear 
from  the  figure  that  under  certain  conditions  no  mode  locking  so¬ 
lutions  are  found.  Indeed,  that  happens  for  a  fixed  excess  gain 


1  -  —  + 


o  “lV 


(2.14) 


parameter  1  -  l/gQ  when  gQic  becomes  too  large,  i.e.  the  Q  of 
the  absorber  becomes  too  small  or  its  saturation  power  becomes 

too  small.  The  saturable  absorber  is  too  "overpowering",  the  pulse 
wants  to  become  too  high  and  too  short,  and  the  laser  medium  cannot 
adjust  itself  to  it. 

An  asymptotic  expression  for  l/wT in  as  a  function  of  g  < 
is  easily  obtained  and  helps  in  the  interpretation  of  Figure  2. 

If  one  assumes  that  go«:  <<  1,  i.e.  if  one  views  the  portion  of 
the  plots  closest  to  the  ordinate,  and  looks  for  solutions  1/wtT,, 
<<  1,  one  may  simplify  (2.14)  to  become 


(2.15) 


Thus  the  lower  branches  of  the  curves  in  Figure  1  are  straight 
lines  with  slopes  proportional  to  [l  -  —  . 

I  V 

Since  we  have  found  in  general  two  solutions  for  each  set  of 
parameters  corresponding  to  different  pulse  widths  and  intensities 
one  would  expect  that  in  general  some  of  these  solutions  would  be 
unstable.  This  question  will  be  investigated  in  another  publica¬ 
tion.  Suffice  it  to  state  here  that  the  branch  above  the  "locus 
of  apices"  in  Figure  1  is  found  to  correspond  to  unstable  solu¬ 
tions,  except  in  the  very  neighborhood  of  the  locus.  There  two 
stable  solutions  are  found  within  a  very  narrow  regime  of  xg^, 
for  any  given  gQ. 


Conclusions : 


The  problem  of  a  fast  saturable  absorber  has  been  solved  in 
closed  form.  The  pulse  shape  is  a  secant  hyperbolic  as  a  function 
of  time.  Incidentally,  the  spectrum  has  also  a  secant  hyperbolic 
shape  because  the  secant  hyperbolic  is  its  own  Fourier  transform. 

The  net  gain  of  the  system  (g  -  1)  .is  negative  before  and 
after  the  pulse.  This  is  a  necessary  requirement  for  the  stability 
of  the  solution  found.  Indeed,  if  the  gain  had  been  found  positive 
in  the  dead-time  between  the  pulses,  any  disturbance  could  grow  in 
this  time  interval  and  thus  destroy  the  sequence  of  periodic  short 
pulses .  Experimentally  observed  two  photon  fluorescence  traces 
produced  by  picosecond  pulses  have  been  compared  against  assumed 
gaussian  and  Lorentzian  pulse  shapes  [1]  .  The  present  analysis 
suggests  that  the  assumption  of  a  hyperbolic  secant  shape  is  more 
plausible . 


